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point  observations  of  single  bathythermographs.  Few  models 
produce  an  accurate  picture  of  the  ocean  environment  that 
can  be  used  for  updatinq  tactics  to  conform  to  a  changing 
situation.  Producing  a  reliable  prediction  of  conditions 
for  a  large  area,  while  using  limited  resources,  is  the 
basic   objective    of   this   paper. 

Satellite  infrared  imaging  of  the  ocean  surface  has  been 
used  effectively  to  map  sea  surface  temperature  patterns. 
Such  sea  surface  temperature  patterns  can  be  used,  along 
with  climatology,  tc  identify  subsurface  thermal  structure 
in  an  ocean  area  according  to  results  of  this  study.  More 
accurate  inputs  can  be  made  to  range  dependent  acoustic 
prediction  models,  thus  improving  the  ASW  environmental 
predictions   available   to   fleet   users. 
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I.       INTRODUCTION 


A.       PREDICTION    OF    OC2AN    THERMAL    STRUCTURE 


Attempts  to  predict  the  vertical  temperature  profile  in 
a  large  body  of  water  can  be  quite  varied.  A  climatoiogical 
atlas  may  be  used  and  a  profile  chosen  which  corresponds  to 
the  location  and  time  of  interest.  The  climatology  may  be 
updated  using  an  interpolation  from  one  climatoiogical 
period  to  another.  A  recent  measurement  in  an  area  can  be 
used  to  modify  the  climatoiogical  prediction,  thus  producing 
a  more  -ime-sensitive  prediction.  Measurements  which  may  be 
used  to  to  modify  climatoiogical  predictions  are  temperature 
structure  in  the  upper  ocean  and  sea  surface  +emperature 
(SST)  . 

The  bathythermograph  (BT)  provides  a  temperature  profile 
at  a  point,  while  satellite  measurements  can  provide  recent 
sea  surface  temperature  patterns  over  an  area  of  interest. 
A  combination  of  climatology,  recently-measured  temperature 
structure  (ET) ,  and  satellite-derived  SST  could  provide  a 
very  reasonable  prediction  of  the  vertical  temperature 
profile   for  a   body   of    water. 

3.       PROBLEM    OF    ACOUSTIC    PREDICTIONS 

Accurate  predictions  cf  transmission  loss  depend  on 
accurate  specifications  of  water  masses  along  an  acoustic 
path.  The      primary    physically-measurable      quantity      which 

affects  sound  speed  (and  thus  transmission  loss)  is  tempera- 
ture. Temperature  structure  of  the  water  is  an  important 
factor  in  acoustic  losses  along  a  transmission  path.  There 
is   no      way    to      gather    temperature      profiles   rapidly      along   a 
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desired   track.  To    measure     temperatures   along      the   tracks 

analyzed  for  the  experiment  discussed  in  this  thesis 
required  a  nine-hour  mission  by  a  dedicated  patrol  aircraft. 
The  operational  fleet  user  does  not  have  the  luxury  of  dedi- 
cated air  assets  to  provide  a  detailed  thermal  description 
of   an   area    of   acoustic  interest. 

Obviously,  the  employment  of  a  dedicated  patrol  aircraft 
loaded  with  a  store  of  expendable  bathythermographs  is  not  a 
realistic  solution  to  the  problem  of  predicting  acoustic 
properties    in   various      locations   of   the    world.  The   use    of 

numerical  models,  relying  on  climatological  files,  and  modi- 
fied by  recent  observations,  provides  a  first  guess  in 
predicting  thermal  structure.  However,  these  models  require 
time  to  produce  results,  and  may  therefore  not  accurately 
present   the  actual    thermal    structure. 

If  a  coupling  exists  between  conditions  at  the  air-sea 
boundary  and  the  depth  of  the  mixed  layer,  satellites  can 
provide  a  means  of  rapidly  surveying  the  sea  surface  in  an 
oceanic     operating        area.  The        satellite      sea        surface 

observations  could  then  be  used  to  estimate  the  vertical 
temperature  structure  throughout  the  operating  area.  If  a 
good  linkage  exists  between  sea  surface  conditions  and  the 
mixed  layer  depth,  the  accuracy  and  timeliness  of  acoustic 
predictions      could        be      improved.  Such      ocean        surface 

conditions  observable  remotely  by  satellite  include  sea 
surface   temperature,    ocean    color,    and    topography. 

C.       EXPERIMENTAL    BASIS    FOR    THIS    THESIS 

In  November  and  December  of  1980,  the  Naval  Postgraduate 
School  conducted  the  Acoustic  Storm  Transfer  and  Response 
Experiment  (ASTREX) .  This  experiment  was  conducted  as  part 
of  the  first  storm  Transfer  and  Response  Experiment  (STREX) , 
a      joint      investigation      by      United      States      and      Canada      to 
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Figure    1.1        Eastern    North    Pacific   with    ASTHEX  Regior 


Defined. 
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determine  the  effects  on  the  air-sea  boundary  of  North 
Pacific  storms,  and  the  converse  [Ref.  1  ].  Figure  1.1  shows 
the  location  of  the  ASTSEX  area  in  the  context  of  the 
eastern      North   Pacific     Ocean.  Three      flight   tracks      were 

examined  within  this  area.  Each  track  covered  a  diagonal 
through  the  central  area  of  the  square  from  the  southeast  to 
the    northwest. 
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II.       CLIMATOLOGY    OF    THE    EASTERN    NORTH    PACIFIC    OCEAN 

A.       OVERVIEW    OF    THE    EASTERN    NORTH    PACIFIC    OCEAN 

Much  of  the  research  conducted  in  the  Eastern  North 
Pacific  Ocean  relies  on  data  collected  at  Ocean  wea-.her 
Station  "P"  (OWS-P) ,  located  at  50N,  145W,  in  the  center  of 
the    Gulf      of    Alaska.  The    research      project   reported      here 

operated  in  an  area  of  the  North  Pacific  •'re  the  east  and 
south  of  OWS-P.  Tabata  [Ref.  2]  concluded  that  the  charac- 
teristics of  OWS-P  were  representative  of  the  waters  in  the 
northeast   Pacific   region 

Ocean  Station  "P"  is  situated  north  of  the  Subarctic 
Current,  well  to  the  northwest  of  the  divergence  which  is 
formed  when  the  West  Wind  Drift  approaches  the  North 
American  continent  [Ref.  3].  The  California  Current  is  the 
southerly  component  of  the  West  Wind  Drift  after  it  divides, 
while  the  northerly  component  forms  "rhe  Alaska  Current  and 
continues   i*s    flew    into    the    Alaska   Gyre   [Ref.    4]. 

The  climat  ological  description  of  •'■he  project  area 
depends  on  a  combination  of  multiple-year  observations  from 
OWS-P  and  observations  within  the  area  itself.  The  region 
from  which  the  data  were  taken  includes  the  subarctic  front 
and    the    source      of    the  California   Current.  The    California 

Current    is    a    permanent  feature    in    the    observation    area. 

The  observation  area  may  be  described  by  three  regions: 
the  Coastal,  or  Transition  Region,  extends  from  the  coast 
out  to  130W;  the  Subarctic  Watermass  extends  seaward  beyond 
the  coastal  zone  to  the  north  of  42N;  and  -he  Central 
Pacific  Watermass  is  south  of  42N.  [Ref.  5  1.  Boundaries 
separating  these  zones  are  dynamic,  and  demarcation  lines 
serve    only    as    initial    positions    for    observing   variability. 
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A  seasonal  change  in  the  location  of  these  boundaries 
can  be  no-ed  from  observations  of  climatological  data. 
These  seasonal  changes  can  be  related,  in  part,  to  weather 
conditions  in  the  Northeast  Pacific.  In  the  winter,  the 
Aleutian  Low  dominates  the  atmospheric  sea  level  pressure 
patterns,  and  in  the  summer  the  region  is  under  the 
influence  of  the  North  Pacific  High.  These  systems  both 
drive  'he  prevailing  westerly  winds,  which  give  rise  to  the 
West  Wind  Drift  in  the  upper  ocean  [Ref.  7]. 

Although  the  boundaries  do  vary,  the  variation  of  the 
boundary  between  the  Subarctic  and  the  Subtropic  watermasses 
is  low.  The  position  of  the  front  was  "relatively 
constant",  between  UON  and  42N,  as  observed  by  Dodimead  et 
al.  [Ref.  8].  Rod  en  [Ref.  9]  found  the  position  of  the 
front  to  vary  from  40 N  to  4  5N,  and  stated  that  this  frontal 
position  was  highly  dependent  on  the  prevailing  wind  fields. 
This  seems  reasonable  for  a  latitudinal  boundary  as  the  wind 
direction  remains  predominately  westerly  from  month  to 
month,  while  wind  speed  in  winter  is  generally  twice  that  of 
summer  [  Ref.  10  ]. 

Dodimeai  e_t  aj..  found  that  the  coastal  zone  was  much 
smaller  in  winter  than  in  summer  [Ref.  11]-  This  could  be 
due  in  part  to  an  increase  in  wind  speed  in  winter  (compared 
to  summer)  along  a  longitud indai  boundary  that  is  terminated 
by  the  coastal  zone  and  North  America. 

B.   WATER  MASS  CHARACTERISTICS 

The  region  used  for  A5TREX  contains  three  water  masses. 
These  are:  1)  Coastal  Water  extending  from  the  coast  to  the 
California  Current  (130W);  2)  Eastern  Subarctic  Pacific 
Water  extendina  southward  to  about  U2N;  and  3)  Subtropic 
Pacific  Water  to  the  south  of  42N  [Ref.  12]. 
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Dodimead  et  al.  reviewed  the  area  in  terms  of  various 
domains.  This  description  shows  three  zones  in  ♦'he  vertical 
salinity  structure:  the  upper  zone,  located  above  a  perma- 
nent halocline,  displays  a  seasonal  variation;  the  permanent 
halocline  (100-200m)  represents  the  second  zone  whose 
vertical  extent  varies  with  latitude;  below  the  halocline  is 
the  lower  zone,  which  begins  at  the  salinity  surface  (S  = 
33.8  g/kg)  ,  and  has  gradually-increasing  salinity  and 
decreasing      temperature      with      depth.  It      displays      less 

seasonal  variation  in  temperature  and  salinity  than  the 
upper   zone.      [Hef.    13]. 

In  the  upper  zone,  the  region  is  influenced  by  both 
Coastal  and  Transitional  Domains.  The  lower  zone  can  be 
partitioned  into  the  Central  Subarctic,  Transitional  and 
California   Undercurrent   Domains  [  Ref .    14], 

The  Tabata  and  Dcdimaad  et  al.  descriptions  vary  in 
approach,  but  results  are  similar  for  the  watermass 
structure  of  the   AST3EX  region. 

1 •      Watermass    Ana  ly^sis 

The  Coastal  Watermass,  referred  to  as  the  Transition 
Region  by  Sverdrup,  Johnson  and  Fleming  (1942),  extends  from 
the  coast  of  California  to  about  130R  Longitude,  where  it  is 
bounded  by  the  Eastern  North  Pacific  Central  Water  (or 
Subtropic   Pacific   Water  according   to   Tully,    1964) . 

Both  the  southward  California  Current  and  the 
inshore        ccuntercurr ent         (Davidson      Current)  make        the 

Transition   Zone      highly   variable.  By    contrast,         bo-h   the 

Subarctic  Pacific  Water  and  the  Subtropic  Pacific  Water  have 
much    less   variability   over    time. 

Predominately  during  spring  and  summer,  upwelling 
occurs   off    the  California   coast.  in    regions   of    upwelling, 

the    sea   surface   temperatures   in    summsr      may   be   lower   than   in 
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Three   Zone   Structure   in   the   Eastern    North 
Pacific   Ocean     [Hef.     15], 
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winter.  The  Davidson  co untercurrent  continues  to  flow 
northward  during  the  upwelling  season,  but  at  depths  greater 
than  200  meters.  During  the  winter,  the  Davidson  Current 
may  extend  to  the  surface  and  flow  northward  to  about  48N. 

An  examination  of  monthly  mean  depths  to  the  *cp  of 
the  thermocline  shows  a  depth  variation  of  15  to  60  meters 
at  the  end  cf  the  upwelling  season  (July-September).  A  vari- 
ation of  45  to  90  meters  in  the  depth  of  the  thermocline  can 
be  noted  at  the  end  of  the  winter  season  (January-February). 
Additionally,  in  all  seasons,  the  mixed  layer  deepens  with 
increasing  distance  from  the  coast.   [ Ref .  16]. 

Close  to  the  coast,  in  the  Transition  Zone,  salinity 
and  temperature  structures  will  be  affected  by  both 
upwelling  and  coastal  run- off.  With  increasing  distance 
from  the  coast,  a  more  stable  structure  develoDs.  West  of 
130W,  salinities  at  the  surface  are  32.5  to  33.0  g/kg, 
increasing  to  about  33.0  to  33.5  g/kg  at  120  meters  depth. 

The  Eastern  Subarctic  Pacific  Wa~ermass  is  charac- 
terized by  three  layers.  The  upper  layer  (surface  to 
approximately  100  meters)  is  seasonal  in  structure.  In 
winter  this  layer  is  well-mixed,  but  in  summer  a  sharp  ther- 
mocline and  a  halocline  form  in  this  layer.  The  second 
layer  is  the  "principal"  halocline  (from  100  to  200  meters) 
which  marks  the  transition  from  the  seasonal  layer  to  the 
"stable"  deep  layer.  In  the  deep  or  lower  layer  (depths 
greater  than  200  meters),  the  salinity  increases  with  depth 
and  the  temperature  decreases.  At  1000  meters  depth,  the 
temperature  is  2.8C  and  salinity  is  34.4  g/kg. 

By  contrast,  the  Subtrcpic  Pacific  Water  has  no 
permanent  halocline.  Rather,  it  shows  a  salinity  decreasing 
with  depth  to  a  minimum  between  200  and  800  meters. 
Temperature  is  the  dominating  factor  in  the  formation  of  the 
seasonal  layer.     The  upper  zone,   which   displays  seasonal 
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temperature  variations,  extends  to  a  depth  of  about  150 
meters,  beyond  which  a  permanent  thermcciine  exists.  This 
seasonal  zone  is  isothermal  to  a  dep^h  of  150  meters  in 
early  spring,  shallowing  to  a  depth  of  40  to  60  meters 
during  the    summer  and    into    fall. 

2.      Domain   Analysis 


By  partitioning  the  area  into  well-defined  domains, 
a  more  detailed  description  of  the  water  masses  can  be  made. 
However,  there  are  few  clear-cut  boundaries  in  either  the 
upper  or  lower  zones,  and  the  halocline  vanishes  in  certain 
situations    [Ref.    17]. 

The  upper  zone  of  the  test  ragion  is  a  combination 
of  Transitional  and  Coastal  Domains.  The  Coastal  Domain  has 
a  thermal  structure  which  is  very  dependent  on  location. 
Such  location-dependent  processes  as  river  run-off,  heating 
and  cooling,  and  wind  cause  the  Coastal  Domain  to  be  highly 
variable   from    place   to   place  [Ref.    18]. 

The  Transitional  Domain  can  be  identified  by  rela- 
tively warm  waters  (T>15C  in  summer  and  T>7c  in  winter),  and 
ralatively  high  salinities  (S>32.2  g/kg  at  the  surface  and 
S>3  3.U   g/kg      at    the      top    of      the    halocline).  In    the      test 

region  however,  the  salinities  may  be  somewhat  lower 
[Ref.  19]  (as  discussed  earlier,  mean  salinities  at  the 
surface  in  the  •test  region  were  32.5  to  33.0  g/kg;  at  the 
top   of   the    halocline    the    salinities    w  =  re    33.0    to    33.5    g/kg) . 

3ased  on  the  domain  description,  the  upper  zone  of 
the  test  region  can  be  characterized  by  a  layer  cf  water 
approximately      one        hundred     meters      deep,  which      varies 

seasonally.  This  upper  zone  shows  an  increase  of  salinity 
with      depth      to    the      halocline      whan      it    exis-.s.  An      even 

stronger  increase  cf  salinity  occurs  through  the  halocline 
(.003   g/kg/m   <   GRAD    ( S)       <       .008    g/kg/m) .         At   the    bottom   cf 
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Figure   2.4        Upper   Zone  Domains    [Ref.    20] 
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the  halocline  the  salinity  approaches  33.8  g/kg.  Also,  a 
temperature  minimum  may  occur  at  the  bottom  of  the  upper 
layer   [Ref.    21]. 

The  upper  zone  responds  to  atmospheric  affects  and 
shows  much  more  variability  than  does  the  lower  zone.  The 
halocline  acts  as  a  "permanent:  boundary"  between  the  upper 
zone  and  the  lower  zone.  Compared  to  the  uppsr  zone,  the 
lower   zone    has   a  stable   salinity   structure.      [Ref.    22]. 

A  second  halocline  forms  in  +he  upper  zone  during 
the  summer.  This  halocline  begins  at  10  *:o  30  meters  depth 
and  is  due  to  the  mixing  effects  of  the  light  summer  winds. 
Below  this  secondary  halocline  (which  can  extend  to  75 
meters),  the  water  is  nearly  isohaline  down  to  -he  principal 
halocline  (100  to  20  0  meters  dep^h  in  the  Subarctic  Water) 
[Ref.    23]. 

The  lower  zone  of  the  test  region  can  be  partitioned 
into  (at  most)  three  domains,  depending  on  the  season. 
These  lower  domains  are  the  Central  Subarctic,  Transitional, 
and   the   California   Undercurrent    [Ref.    24]. 

The  Central  Subarctic  Domain  is  typical  of  the  lower 
layer  water  in  the  northwest  cornsr  of  the  test  region. 
This  is  also  the  water  type  associated  with  OWS-P.  Dodimead 
et  al.  identify  this  domain  as  water  underlying  a 
well-defined  halocline.  At  OWS-P,  the  halocline  is  defined 
by  a  salinity  gradient  of  0.91  g/kg  per  meter.  The  salinity 
increases  with  depth  very  slowly  in  the  lower  layer 
[Ref.    25]. 

The  Transitional  Domain  underlies  its  surface 
counterpart  throughout  the  southern  half  of  the  test  region 
and  also  extends  through  the  northeast  quarter  of  the 
region.  Near  the  coast,  both  this  domain  and  the  California 
Undercurrent  can  exist.  The  Transitional  Domain  is  strongly 
influenced    by      the    major   geostrophic      currents   in      the    area. 
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Figure    2.6        Lower    Zone   Domains     [Ref.    27] 


28 


It  is  a  narrow  band  underlying  the  West  Wind  Drift  and 
broadens  as  this  current  separates  tc  form  the  Alaskan  and 
California      currents.  In      the        Transitional      Domain      the 

halocline  is  less  pronounced  than  in  the  Central  Subarctic 
Domain.  The  Transitional  Domain  bridges  the  gap  between  the 
strongly-layered  Subarctic  Water  and  the  Subtropical  Water 
which   has   no   distinct    halocline  (Ref.    28], 

The  third  lower  domain  present  in  the  test  ar<=a  is 
the  California  Undercurrent  Domain.  As  the  name  implies, 
this  domain  is  influenced  principally  by  the  subsurface 
countercurrent  which  runs  northward  near  the  coast.  The 
countercurrent  originates  in  the  warm  wat<=r  of  the  Eastern 
Central  Pacific  Water.  Dodimead  e_t  aj^.  characterizes 
California  Undercurrent  water  by  temperatures  greater  than 
6C  on  ■khe  34  g/kg  surface  of  salinity,  as  opposed  to  the 
Transitional  Domain  defined  at  a  salinity  of  33.8  g/kg  with 
temperatures   less    than  6.5C    [Ref.    29]. 

3  .      Temperature   S  truct  ure 

The  salinity  structure  of  the  eastern  North  Pacific 
provides  a  relatively  stable  basis  for  defining  water  masses 
in  the  area.  However,  characterization  cf  water  masses 
requires  a  salinity-temperature  relationship  as  well.  The 
true  variability  of  acoustic  conditions  in  the  eastern  North 
Pacific  must  come  from  periodic  and/or  random  variations  in 
temperature  which  are  superimposed  on  the  basic  density 
structure      determined    by      salinity.  Of   course,         salinity 

variations  may  be  of  importance  in  the  coastal  regions  too, 
where  siginificant  intrusions  of  fresh  water  or  more  saline 
upwelled   water   may    occur. 

At  OWS-P,  the  temperature  of  the  upper  zone  has  an 
annual  periodic  variation,  reaching  a  maximum  in  the  late 
summer,      and  a   minimum   in      late   winter.         The    maximum   annual 
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temperature  variation  notsd  between  1956  and  1962  was  about 
9C.  Variations  of  temperature  in  the  entire  upper  zone 
followed  the  same  cycle  as  sea  surface  temperature, 
appearing  to  be  in  phase  with  the  periodic  sea  surface 
temperature  variation  [Ref.  30]. 

As  with  salinity,  the  temperature  variations  are 
greatest  in  the  upper  zone.  The  temperature  range  decreases 
from  a  maximum  of  9C  in  the  upper  53  meters  to  less  than  3C 
variation  at  125  meters  depth.  Below  the  halocline,  the 
small  temperature  variation  with  depth  (gradient)  should  be 
reasonably  well  predicted  from  monthly  climatological  data. 

In  the  upper  zone,  a  seasonal  thermocline  begins  to 
develop  in  spring  with  the  advent  of  the  heating  season. 
The  isothermal  March  layer  becomes  a  two-layer  medium 
through  April  when  a  seasonal  thermocline  begins  to  develop 
at  abcu*  30  meters.  This  seasonal  thermocline  reaches  a 
maximum  depth  of  about  50  meters,  depending  on  the  strength 
of  the  prevailing  winds,  but  can  be  expected  to  occur 
between  30  and  50  meters  during  the  spring  and  summer.  The 
thermocline  shews  the  strongest  gradient  in  late  summer 
(August-September),  when  it  is  at  its  shallowest  depth  (30 
meters)  [Ref.  31  J. 

With  the  advent  of  autumn,  the  cooling  season 
begins.  The  seasonal  thermocline  deepens  throughout  the 
cooling  season  due  to  the  effects  of  colder  surface  water 
and  increased  mixing  caused  by  the  strong  winter  winds. 

The  seasonal  thermocline  forms,  coincident  with  the 
formation  of  the  secondary  or  seasoral  halocline.  This 
thermocline  is  very  stable,  and  hence,  the  structure  of  the 
upper  zone  is  quite  predictable  [Ref.  32].  The  seasonal 
thermocline  and  permanent  halocline  locations  can  be  used  as 
inputs  to  predict  the  vertical  structure  of  the  upper  100 
meter  layer  during  the  heating   and  cooling  seasons.    Sea 
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surface  temperatures  may  provide  input  regarding  the 
strength  of  the  seasonal  thermocline  (gradient)  since  both 
the  layer  thickness  and  the  temperature  at  the  bottom  of  the 
thermocline  remain  relatively  constant. 

Characteristics  of  the  Subarctic  Water  are  found  in 
the  waters  which  comprise  the  Transitional  Domain.  Waters 
in  the  Coastal  Domain,  on  the  other  hand,  will  vary  from  the 
Subarctic  Water  structure.  Coastal  Waters  remain  very  near 
the  coast  in  winter,  but  extend  further  seaward  in  the 
summer.  The  primary  effect  of  mingling  Coastal  Waters  with 
Transition  Water  will  be  on  salinity  structure.  The  Coastal 
Water  rends  to  be  less  saline  than  Transition  Water.  A 
winter  section  of  temperature  and  salinity  taken  through  the 
test  region  in  1959  indicated  a  narrow  zone  of  coastal 
influence  (*-.o  about  130W)  ,  which  showed  as  an  increase  in 
temperature  and  decrease  in  salinity  from  the  Transitional 
Domain.  A  summer  section  at  the  same  latitude  showed  a 
strong  coastal  influence  out  to  130W;  the  temperatures  were 
much  warmer  and  salinities  much  lower  in  the  Coastal  Domain 
than  in  the  Transitional  Domain. 


C.   METEOROLOGY  OF  THE  EAST EPN  NORTH  PACIFIC  OCEAN 

The  discussion  of  the  meteorology  of  -"-he  ASTREX  region 
and  the  northeast  Pacific  Ocean  will  be  limited  to  factors 
which  affect  spaceborne  instruments,  rather  than  those 
factors  which  influence  the  circulation  and  heat  transfer  in 
the  area. 

The  most  obvious  path  interference  encountered  by  a 
spaceborne  instrument  is  atmospheric  moisture  in  the  form  of 
clouds.  During  ASTREX,  clouds  obscured  most  of  the  satel- 
lite field  of  view  on  every  measurement  day.  Since  the 
purpose  of  the  experiment  involved  studying  the  effects  of 
North  Pacific  storms,  cloudiness  was  an  unavoidable  obstacle 
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for  the  secondary  mission  of  comparing  satellite  data  to  In 
Situ  observations.  Winter  cloudiness  in  the  North  Pacific 
is  a  major  problem  for  visual  and  infrared  imaging  of  the 
ocean  surface. 

The  meteorological  phenomena  which  directly  affect 
space-based  ocean  viewing  instruments  are  aerosols  in  the 
form  of  moisture  and  salts.  Moisture,  in  its  obvious  forms 
of  clouds  and  rain,  severely  degrades  passive  observation 
instruments  which  rely  on  visual  or  infrared  energy. 
Passive  microwave  radiometers  suffer  degradation  when  the 
viewing  area  includes  sevsre  rains. 

The  Navy  Marine  Climatic  Atlas  of  the  World  [  Ref .  33] 
was  consulted  to  determine  climatic  conditions  in  the 
Eastern  North  Pacific  in  November  and  December.  For  both 
months,  cloudiness  greater  than  4/8  can  be  expected  about 
four  days  in  five.  This  represents  a  severe  degradation  for 
visual  and  infrared  imaging  cf  the  North  Pacific  waters  in 
November  and  December. 

Precipitation  reported  in  the  Atlas  increases  from  19 
percent  in  November  to  about  2  4  percent  in  December.  In 
other  words,  rain  occurs  about  one  day  in  fiva  in  November 
and  one  day  in  four  in  December.  Rain  constitutes  an  inter- 
ference effect  for  passive  microwave  systems  and  total 
degradation  for  visual  and  infrared  systems. 

Relative  humidity  was  reported  as  greater  than  90 
percent  one  day  in  five  in  November  and  one  day  in  four  in 
December.  This  follows  the  precipitation  variation.  Maul 
and  Sidran  [Ref.  34 ]  reported  on  effects  of  water  vapor  at 
45N  and  60N.  From  their  report,  water  vapor  content  in  the 
atmosphere  over  the  Eastern  North  Pacific  (between  4Q-50N) 
can  cause  a  temperature  error  of  the  order  of  2  degrees 
Kelvin  in  winter  and  6  degrees  Kelvin  in  summer.  These 
errors   represent    the   difference    between   uncorrected 
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satellite    sea      surface  temperature      observations    and      actual 
sea    surface  temperatures. 
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III.       SATELLITES    IN    OCEANOGRAPHY 

A.       INSTRUMENTS    FOR    OBSERVING    THE   OCEANS 
1 •      Visual   and   Infrared    Imaging 

Several  series  of  satellites  have  both  visual  and 
infrared  capabilities.  Among  these  satellites  are  the  GOES 
series,      NIMBUS      series,      and   NOAA   series.  Stewart    (1981) 

presents  an  excellent  summary  of  satellite  instruments  used 
for  oceanography.  He  also  lists  the  satellites  which  carry 
oceanographic  instruments  and  the  specific  instruments  each 
satellite   carries   [Ref.    35]. 

The        Geostationary  Observational        Environmental 

Satellite  (GOES)  system  is  composed  of  two  satellites  in 
geosynchronous   earth      orbits  over   the    equator.  The   satel- 

lites, designated  GOES-East  (subpoint  at  75  degrees  west 
Longitude)  and  GOES-West  (subpoint  at  135  degrees  west 
Longitude) ,  scan  the  earth,  with  overlapping  coverage  over 
the  Americas,  from  Africa  to  the  western  Pacific  Ocean 
[Ref.    36]. 

The  GOES  satellites,  using  the  visible  and  infrared 
spin  scan  radiometer  (VISSR),  can  scan  the  full  disk  of  the 
earth  in  18.2  minutes  and  generate  an  image  every  thirty 
minutes  [Ref.  37].  Figure  3. 1  is  a  visual  image  produced  by 
the  GOES  VISSR.  Table  I  lists  the  spectral  channels  and 
resolutions  for  the  VISSR.  Figure  3.2  presents  the  atmos- 
pheric transmission  windows  for  the  visible  and  infrared 
portions  of  the   electromagnetic    spectrum. 

GOES  can  provide  near-real-time  coverage  of  oceanic 
areas  in  the  Atlantic  and  Pacific  Oceans.  GOES  imagery  was 
utilized      during    the      STREX    experiment      as    a      meteorological 


3a 


r ABLE    I 

Visible  and   Infrared   Spin  Scan   Radiometer    (VISSR) 

Channelization. 


Channel        Spatial   Resolution  Wavelength 

{at   subpoint)  (micrometers) 

Visible  1    km  0.55    -    0.75 

Infrared  8    km  10.5-12.6 

Note:  .  .    . 

1.      Fett   et   aj..    list   spatial   resolution    for   the   visible 

channel  as   0.5    nautical    miles,    and   5    nautical    mil?s 

for   the   infrared   channel. 


aid.  GOES  images  of  sea  surface  temperature  are  useful  for 
studying  oceanic  circulation,  but  the  spatial  resolution  is 
not  as  good  (8  km)  as  ths  resolution  from  satellites  in  low 
ear^h  orbits,  such  as  NOAA-6  (1.1  km  resolution  at  infrared 
wavelengths)  . 

The  Nimbus-7  spacecraft  carries  the  Coastal  Zone 
Color  Scanner  (CZCS) ,  which  is  a  visual  radiation  instrument 
centered  in  specific  wavelength  bands  to  measure 
concentrations  of  phytoplankton  in  the  ocean.  Table  II 
shows    the    wavelengths    for   indicated   CZCS    channels. 

The  Advanced  Very  High  Resolution  Radiometer  (AVHR3) 
is  a  five-channel  instrument  used  on  the  NOAA-6  and  NOAA-7 
spacecraft.  The  AVHRR  on  the  NOAA-6  has  one  visual  channel 
and  three  infrared  channels  (the  fifth  channel  is  a  repeat 
of  channel  four)  .  NOAA-7  utilizes  all  five  channels.  The 
AVHRR  has  a  ground  resolution  of  1.1  km  in  ail  channels.  In 
the  daytime,  the  first  two  channels  (visible  and  far  red  to 
near  IR)  can  be  used  to  filter  cloud-contaminated  infrared 
measurements.  Channel  three  (3.5-3.93  um)  can  only  be  used 
at  night  because  of  severe  contamination  by  reflected  solar 
energy.  Channel  four  (10.5-11.5am)  can  be  used  day  or  night 
for    measuring   sea    surface   temperature. 
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Figure    3.1        GOES    Visible  Image    (Full   Disk)     [Ref.    38], 
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Figure  3.2 


Atmospheric  Tr ansmittance  (Visible  and  Infrared 
Spectra)  [Ref.  39], 
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TABLE  II 
Coastal  Zone  Color  Scanner  (CZCS)  Channelization, 

Channel      Wavelength  (microne  ters) 

1  0.433  -  0-453  (blue) 

2  0.510  -  0.530  (green) 

3  0.540  -  0.560  (green 

4  0.660  -  0.680  (red) 

5  0.700  -  0.800  (far  red/near  IR) 

6  10.5-12.5    (far  IR) 

Note: 

1.  Channel  6  failed  in  197  9. 

2.  Resolution  is  approximately  1  km  for  all  channels. 


Table  III  is  adapted  from  both  Reference  36  and  37. 
It  lists  the  AVHRR  spectral  channels  fcr  the  AVHRR/1  and 
AVHRR/2. 


TABLE  III 

Advanced  Very   High   Resolution   Radiometer    (AVHRR) 

Chan  nslization 


Channel  Wavelength    (micrometers) 

1  0.58   -    0.68         (visible) 

2  0.725    -    1.10      (far   red-nsar    IR) 

3  3.55    -    3.93  middle    IR) 


4  10.5    -    11.5  far    IR) 

5  1 1. 5   -    12.5         (far   IR) 

Note: 

1.  Spatial   resolution   at   satellite    ?ubpoint    is 
0.5    nautical   miles    (1.1    km). 

2.  Channel   4    wavelenath    is    10.3    -    11.3    urn    for   all 
AVHRR/2    instruments. 

3.  Channel   5   added   to   further   enhance    sea    surface 
measurements   in   "he   tropics. 
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The  noise  equivalent  temperature  (NETD)  of  the  AVHRR 
is  less  than  0.12  degrees  Kelvin  for  a  300  degree  Kelvin 
scene  (Ref.  40].  The  temperature  sensitivity  for  the  chan- 
nels is  approximately  0. 5 -degrees-Kelvin  when  viewing  a 
300-degree-Kelvin  scene.  Since  the  noise  equivalent  temper- 
ature of  the  system  is  less  than  the  temperature  sensitivity 
of  the  channels,  identification  of  thermal  gradients  on  thQ 
order   of   0 . 5-degrees-Kelvin    per    kilometer   is    feasible. 

In  Chapter  two,  the  meteorological  phenomena  which 
interfere  with  visible  and  infrared  imagery  of  the  ocean 
surface  were  discussed.  The  following  discussion  deals  with 
the  AVHRR  and  techniques  used  to  compensate  for  atmospheric 
effects.  The  AVHRR  measures  a  "brightness  temperature"  of 
the  sea  surface,  and  the  actual  temperature  can  be  computed 
using  physical  radiation  laws.  The  sea  surface  temperature 
is  calculated  from  the  brightness  temperature  by  correcting 
for      attenuation      and      emission      of      the      atmosphere.  The 

radiometer  is  calibrated  on  board  the  spacecraft  by  having 
the  instrument  view  deep  space  (approximately 

2-degree-Kelvin  black  body),  and  the  radiometer  housing 
(monitored  by  thermistors)  .  The  multi-channel  simultaneous 
imaging  of  the  sea  surface  forms  the  basis  for  the  atmos- 
pheric  corrections    used    for    the    AVHRR    measurements. 

Deschamps  and  Phulpin  [Ref.  41]  proposed  a 
multispec- ral  correction  method  which  is  the  basis  for  the 
AVHRR  measurement  accuracy  of  sea  surface  temperature.  The 
simultaneous        solution        of  three         different        "window" 

measurements  removes  many  non-linearities  caused  by 
atmospheric  effects  and  approximations  used  to  calculate 
radiances.  The  cost  of  improving  atmospheric  corrections  is 
an  increase  in  the  noise  equivalent  temperature  of  the 
system. 
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The  AVHRR/1  (NOAA-6)  uses  two  bands  for  analysis  and 
atmospheric  corrections.  The    analysis    bands      are  channels 

three   and   four.  Each   channel   has   a      different    response    to 

water  vapor  in  the  atmosphere.  The  variation  of  response  of 
the  two  channels  makes  the  corrections  mors  reliable. 
However,  channel  three  operates  very  near  the  five 
micrometer  wavelength  where  solar  energy  is  about  equal  to 
the  emitted  energy  of  a  3 OO-degree-Kelvin  black  body  (sea 
surface).  Therefore,  channel  three  is  a  night  channel. 
Channel    four   can   be   used   day  or    night. 

Leitao  et  al.  [Ref.  42]  discuss  the  use  of  satel- 
lite infrared  imagery  to  locate  the  northern  surface 
boundary  of  the  Gulf  Stream.  They  conclude  that  the  maximum 
sea  surface  temperature  gradient  is  a  good  approximation  to 
the    location   of   the   northern   boundary. 

A  problem  with  determining  sea  surface  temperature 
gradients  from  space  is  the  variability  of  the  atmosphere. 
Maul  [Ref.  43  ]  discusses  the  effects  of  the  atmosphere  on 
ocean  surface  thermal  gradients.  In  a  moist  atmosphere,  a 
thermal  gradient  can  be  reduced  to  less  than  half  of  its 
actual  value   by   the   atmosphere. 

2 •      Passive    Micro  wave   R adio metr v 

The  latest  instrument  used  for  passive  microwave 
observations  of  the  earth  is  the  Scanning  Multi-frequency 
(channel)  Microwave  Radiometer  (SMMR) .  This  instrument  is 
operational  on  NIMBUS-7.  sea  surface  temperatures  are 
calculated  from  brightness  temperatures  measured  at  6.6 
gigahertz  (4.54  cm  wavelength).  At  this  frequency,  spatial 
resolution      is  generally      considered    to      be      150    km.  Test 

results  of  the  HIM5as-7  system  indicated  that  resolution 
cells    for    the   6.6    GHz    channel    were    157x156    km   [Ref.    44]. 
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Lipes  [Ref.  45]  referred  to  the  detrimental  effects 
of  rain  and  sunglint  (reflection  of  the  sun  off  of  the  water 
surface)  on  the  SEASAT  SMMR  (a  similar  instrument  is  on 
NIMBDS-7).  Rain  and  sunglint  caused  the  SMMR  derived  sea 
surface  temperature  to  be  much  higher  than  the  actual  sea 
surface  temperature.  Areas  without  sunglint  or  rain  were 
more  accurately  measured  by  tha  SEASAT  SMMR. 

In  addition  to  problems  with  rain  and  sunglint,  the 
SMMR  is  also  susceptible  to  radio  frequency  interference 
(RFI)  ,  and  severe  sidelobe  contamination  within  600  km  of 
land  [Ref.  46  ].  Improvements  in  the  instrument  design  and 
ir.  the  ground-based  processing  design  would  reduce  the  side- 
lobe  contamination  problem  for  the  SMMR  [Ref.  47], 

3  .   Sa  t  ellite  Radar  Altimeter 

Satellite  altimetry  presents  an  active  approach  to 
locating  water  mass  boundaries  and  other  density  related 
anomalies  in  the  oceans.  The  satellite  measures  height 
above  the  surface  of  the  sar^h.  These  daxaf  combined  with 
the  orbital  parameters  of  the  spacecraft,  form  the  basis  for 
measuring  the  height  of  the  sea  surface  along  the 
sub-satellite  track.  I  he  radar  altimeter  maps  the 
topography  of  the  ocean  along  the  satellite  ground  track, 
since  the  instrument  is  fixed  at  nadir  and  does  not  scan 
across  track.  The  ocean  topography  is  the  difference 
between  the  sea  surface  and  the  geoid.  Surface  gees-trophic 
currents  can  be  derived  from  the  ocean  topography  generated 
by  a  radar  altimeter. 

Altimeter  height  measurement  uncertainties  exist 
due  to  atmospheric  propagation  and  ionospheric  effects, 
local  deviations  of  the  ocean  surface  topography  from  the 
geoid,  and  inaccuracies  in  the  satellite  orbit.  Density 
gradients  associated  with  geostrophio  currents  account  for  a 
small  part  of  the  topographic  variations  [Ref.  48]. 
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The  SEASAT  altimeter,  accurate  to  10  centime-  ers  in 
absolute  altitude  measurene nt ,  was  an  excellent  instrument 
for  measuring  the  location  of  surface  geostrophic  currents. 
The  GEOS-3  altimeter,  accurate  to  50  centimeters  in  absolute 
altitude,  was  capable  of  locating  the  slope  associated  with 
the  Gulf  Stream  by  using  an  averaging  process  [Ref.  49]. 

The  altimeter  is  a  promising  instrument  for  the 
future.  It  is  capable  of  detecting  western  boundary 
currents  such  as  the  Kuroshio,  but  the  error  of  the  height 
measurement  is  still  larger  than  the  signal  generated  by  a 
weak  current.  A  western  boundary  current  will  cause  the  sea 
surface  to  rise  one  meter  over  100  kilometers.  A  10  centi- 
meter accuracy  of  height  measurement  is  required  to  measure 
such  a  slope  of  the  sea  surface  accurately.  When  the  slope 
of  the  sea  surface  is  an  order  of  magnitude  less  than 
1m/100km,  mere  accurate  height  measurements  are  required. 

B.   OTHER  INSTRUMENTS  FOR  OBSERVING  THE  OCEANS 

This  chapter  has  been  directed  toward  those  instruments 
which  may  be  of  importance  to  the  problem  of  predicting  the 
mixed  layer  depth  from  sea  surface  temperature.  The  problem 
is  much  more  complex  than  cause  and  effect  between  sea 
surface  temperature  and  mixed  layer  depth.  This  report  is 
limiting  the  scope  of  the  investigation  to  ascertaining  a 
relationship  between  sea  surface  temperature,  thermocline 
gradient,  and  mixed  layer  depth,  if  indeed  one  exists. 

The  problem  can  be  expanded  to  include  the  effects  of 
wind  mixing  and  heat  budget.  The  space  technology  needed  to 
provide  observations  of  wind  speed,  wave  height,  and  temper- 
ature has  been  demonstrated  on  SEASAT.  The  SMMR  can  provide 
information  on  the  wind  speed,  in  addition  to  sea  surface 
temperatures.  Besides  surface  topography,  the  altimeter  en 
SEASAT  could  provide  information  on  significant  wave  height. 


42 


Additional  instruments  on  SEASAT  provided  wind  velocity 
(Scatterometer)  ,  and  surface  wave  direction  and  wavelength 
(SAR) .  The  Scatterometer  (SCAI)  and  the  Synthetic  Aperature 
Radar  (SAR)  are  active  microwave  imagers  and  not  subject  to 
many   of    the   problems    encountered   by   passive   systems. 
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IV.       EXPERIMENTAL    DESIGN 

A.       OBJECTIVE 

The  objective  of  this  investigation  was  to  determine  if 
a  relationship  between  sea  surface  temperature,  mixed  layer 
depth,  and  thermocline  gradient  ccald  be  established.  By 
reviewing  the  climatology  of  the  area,  and  including  the  sea 
surface  temperature  patterns  and  sub-surface  structure,  an 
empirical  approach  to  forecasting  the  extent  of  the  mixed 
layer  from  surface  observations  and  limited  bathythermo- 
graphs  was    to   be    developed. 

Legeckis  and  Gordon  [Ref.    50]  found  that 

sat  el  lite-  derived  sea  surface  temperature  fronts 

corresponded  to  deep-reaching  subsurface  temperature 
structure.  They   also     concluded   that      the   observed      mixed 

layer  depth  correlated  with  the  surface  temperature  pa-terns 
observed   for   the   Brazil   Current    and   two    warm   core    eddies. 

Blackstone  and  Whritner  [Ref.  51  ]  studied  the  persis- 
tence of  observed  sea  surface  temperature  patterns  to  a 
depth  of  1000  feet  off  the  coast:  of  southern  California. 
Air-dropped        Expendable      Bathythermographs         (AXBT)  were 

employed  to  examine  regions  which  displayed  a  strong  sea 
surface  temperature  gradient  in  Defense  Meteorological 
Satellite  imagery.  The  AXBT's  provided  observations  of  the 
vertical  temperature  structure  associated  with  the 
observable  sea  surface  temperature  patterns.  From  analysis 
of  temperature  patterns  at  various  depths,  whritner  and 
Blackstone  observed  that  the  temperature  gradient  associated 
with  sea  surface  temperature  patterns  was  still  recognizable 
at    10  00    feet. 
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Leitao    et      al.         [Ref.     52]      found    that      the  maximum      sea 

surface   temperature      gradient    (as   observed    from  space)       was 

useful  for  approximating  the  northern  surface  boundary  of 
the    Gulf   Stream. 

B.       DATA    ACQUISITION 

In  November  and  December  of  1980,  the  Naval  Postgraduate 
School  conducted  the  Acoustic  Storm  Transfer  and  Response 
Experiment  (discussed  in  chapter  1)  .  As  part  of  this  exper- 
iment, Air-Dropped  Expendable  Bathythermographs  (AXBT)  were 
placed  along  the  same  track  on  six  days.  There  were  three 
flights      in      November,        covering     the      center      track.  In 

December,  another  three  flights  were  conducted.  The  center 
track  was  again  covered  three  times,  but  on  the  last  two 
flights,  tracks  offset  from  the  center  also  were  flown.  The 
offset  tracks  necessitated  a  reduction  in  the  number  of 
stations  which  could  be  taken  along  either  the  center  track 
or   the   offset   track   [Ref.    53]. 

Satellite  observations  of  sea  surface  temperature  were 
collected  from  infrared  imaging  by  the  Advanced  Very  High 
Resolution         Radiometer      (AVHRR)  en         beard      the        NOAA-6 

spacecraft.  AXBT's  were  placed  along  a  1000-km  track  (at 
55-km  intervals).  The  data  from  the  AXBT's  were  digitized 
in  five-meter  increments  of  depth  (see  Figures  4.1  and  4.2), 
and  the  printed  output  was  used  to  construct  the  vertical 
cross-sections   of   temperature  used   in    this    thesis. 

In  addition  to  -temperature  soundings  collected  from  the 
AXBT's,  soundings  were  made  from  a  research  vessel  in  the 
operating  area  using  a  Neil  Brown  conductivity  temperature 
depth    (CTD)    microprof iler  [Ref.    56]. 
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Figure  4.1 


?Bef ^itf.  T6mPeratur9s   at    Five   Meter   Increments 
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Figure    U.  2        Digitized    Depths  of    Isotherms    [Ref.    55]„ 
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Figure   4.3 


FNOC    Salinity    Profile    Along    the   Flight   Track 
for    13   November. 


48 


C.       DATA    ANALYSIS 
1 .      Da t a    Plots 

As  a  means  of  beginning  an  analysis  of  the  possible 
relationship  between  sea  surface  temperature  patterns  and 
corresponding  patterns  of  mixed  layer  depth,  vertical 
cross-sections  of  temperature  for  each  flight  track  in 
ASTREX  were  constructed.  In  addition  to  these  transects, 
sea  surface  temperature  was  plotted  by  position  along  the 
flight  track  as  well  as  mixed  layer  depth.  A  correlation 
was  no+ed  visually  between  track  position  and  both  sea 
surface   temperature   and   mixed  layer   depth. 

The  measurements  taken  with  the  CTD  were  in  the 
northwest  ccrner  of  the  ASTREX  region.  A  Fleet  Numerical 
Oceanography  Center  salinity  profile  for  the  track  (Figure 
4.3)  shews  very  little  variation  of  salinity  profiles 
between      stations      2      and    13.  Since      the      variability     of 

salinity  in  the  Northeast  Pacific  is  low,  the  soundings 
which  were  taken  with  the  CTD  were  assumed  to  be  reasonably 
representative   of   the    area. 

2«      Statistical   Procedure 

Sea  surface  temperatures  and  mixed  layer  depths  were 
taken  from  Kilonski's  [Ref.  57]  digitized  depths  of  the 
isotherms.  The  mixed  layer  depth  at  a  location  as  reported 
by  Kilcnski  sometimes  shewed  a  variation  of  one  meter 
between  the  two  data  formats  (digitized  temperatures  at 
five-meter      increments        vs.  digitized        depths      of        the 

isotherms)  . 

To    examine    variations  with    position    along    the    track, 
the    outermost      station   on      each    flight      was   assigned      as   the 
origin.         All   other   stations    were      listed   by    (nautical    mile) 
distance   from   the    origin.       For    example,    station    nine    was    the 
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outermost      station    on      19   November.  It      was   assigned      the 

origin  position  (0.0)  for  that  flight  track.  The  adjacent 
station  (seventeen)  is  thirty  nautical  miles  from  station 
nine.      Therefore,    its    position    is   listed   as   thirty    (30). 

In  addition  to  position,  sea  surface  temperature, 
and  mixed  layer  depth,  the  gradient  at  the  thermocline  was 
considered  a  necessary  factor  to  include  in  the  lis-  of 
variables.  Four  depth  intervals  were  selected  for  computing 
gradients   in    the   thermocline.  The    variations   of   gradients 

with  track  position  are  plotted  in  Appendix  A  (November 
tracks)    and   Appendix    B    (December   tracks)  . 

From  the  vertical  cross-sections,  an  initially  small 
interval  was  considered  necessary  for  computing  the 
gradient.  The  smallest  interval  available  from  Kilonski's 
data  was  five  meters.  Since  the  mixed  layer  depths  did  not 
fall  on  multiples  of  five  meters,  the  first  level  above  the 
reported  mixed  layer  depth  was  used  as  a  reference  level  for 
computing  temperature  changes  at  the  top  of  the  thermocline. 
For  example,  if  the  mixed  layer  Slepth  was  67  meters,  the 
temperature  at  65  meters  iepth  was  used  as  the  temperature 
at  the  top  of  the  thermocline.  To  smooth  gradient  varia- 
tions caused  by  this  method  of  defining  the  temperature  at 
the  top  of  the  thermocline,  gradients  were  also  computed 
over    intervals   of   fifteen,     fifty,    and    one-hundred    meters. 

The  computed  temperature  changes  are  listed  in 
various      tables   associated      with      each      flight    track.  The 

normalized  values  of  gradient  computed  from  each  interval 
are    listed    in   Appendix  c. 

The  data  analysis  was  done  using  the  interactive 
statistical      computing      system      MINITAB      [ Ref .    58].  Basic 

statistics  (means  and  standard  deviations),  correlation,  and 
linear  regression  analysis  were  performed  for  each  flight 
track.  MINITAB      correlation      of        data      resulted      in      the 
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computation  of  the  correlation  coefficient  (z)  as  discussed 
in   Miller   and   Freund   [Ref.     59]. 

Three  models  of  mixed  layer  depth  were  developed  for 
each      flight    track.  These   models      ware      simply    the      least 

squares  regression  cf  mixed  layer  depth  as  a  function  of  sea 
surface  temperature,  position,  and  gradient  in  the  thermo- 
cline.  If  the  correlation  between  sea  surface  temperature 
and  mixed  layer  depth  was  weak,  -hen  the  model  of  mixed 
layer  depth  as  a  function  of  sea  surface  temperature  was 
ignored. 

Regression  analysis  in  MINITAB  yields  a  curve  and  a 
coefficient  of  determination  as  a  measure  cf  the  goodness  of 
the  curve-fit  (the  coefficient  of  determination  is  the 
square  cf  the  correlation  coefficient).  To  avoid  confusion, 
each  of  the  measures  of  correlation  between  various  combina- 
tions  of   parameters   will   be    identified    whenever    used. 

D.       OBSERVATIONS    BY    FLIGHT 
1 .      Flight  Track    On a 

The  vertical  temperature  cross-section  for  flight 
track  one  is  shown  in  Figure  4.4.  Ten  stations  were  used 
to  construct,  this  profile.  Isotherms  tended  to  be  deeper 
nearer   the    coast   than   out   at   sea. 

A  warm  water  pocket  exists  between  stations  four  and 
two.  This  warm  pocket  persisted  through  the  vertical  extent 
of   the  soundings   between   these    stations. 

Figure  4.5  shows  both  the  sea  surface  temperature 
and  mixed  layer  depth  along  the  flight  track.  A  possible 
correlation  between  sea  surface  temperature  and  mixed  layer 
depth  along  the  flight  track  is  suggested  by  the  variations 
of   both    with    position    along    the   track     (r    =    0.33). 
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Figure   4. a 


Vertical   Temperature   Cross-section 
(15    November)  . 
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Figure   4.5 


Mixed   Layer   Depth  and   Ssa    Surface   Temperature 
(15   November)  . 
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For  flight  one,  a  strong  negative  correlation  (r  = 
-0.818)  exists  for  the  gradient  computed  from  the  upper  five 
meters  of  the  thermocline.  This  suggests  that  the  initial 
gradient  at  the  thermocline  weakens  in  the  transition  region 
(coastal   waters).  The   weakening  of   this      initial   gradient 

could  also  be  due  to  the  presence  of  a  warm  meander  or  eddy 
as   suggested   by   the  vertical   temperature   cross-section. 

Table   IV    lists,    by    station: 

1.  position   relative   to    the    most-seaward   buoy; 

2.  mixed   layer    depth; 

3.  sea  surface   temperature; 

4.  temperature    change      in   the      upper    five      meters    of      the 
thermocline    (dT/5)  ; 

5.  temperature      change   in     the      upper      15    meters      of      the 
thermocline    (dT/15)  ; 

6.  temperature      change  in     the      upper      50   meters      of      the 
thermocline    (dT/50) ;    and 

7.  temperature      change   in     the   upper      100      meters    of      the 
thermocline    (dT/100)  . 

Table  IV  is  the  matrix  which  was  manipulated  to  examine 
correlations  between  various  elements  in  the  problem,  and 
was  a  source  of  values  in  regression  analysis  for  model 
equations. 

Attempts  to  fit  several  equations  for  mixed  layer 
depth  were  made  using  various  measurable  quantities.  The 
results   are   listed    in    Appendix   D.  Table   V    lists    -"-he    model 

variables  used  and  provides  the  computed  coefficient  of 
determination    for   each   model.  The    most    striking   relation- 

ship occurred  for  a  fit  of  mixed  layer  depth  as  a  function 
of  track  position,  sea  surface  temperature,  and  gradient 
from  100  meters.  The  coefficient  Df  determination  for  the 
fit  was  0.38.  The  correlation  (r)  between  MLD  and  the  func- 
tion     of      the      three      variables      was      0.94.  Position      and 
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T  ABLE    IV 
MLD,    SST   and   Gradients   at   the  Thermooline,    15    November    1980, 


luoy 

Position 

MLD 

SST 

dT/5m 

dT/15m 

dT/50m 

dT/100m 

2 

54  0 

32 

13.81 

1.15 

2.86 

3.90 

4.65 

3 

48  0 

69 

15.88 

1.41 

2.69 

4.62 

5.92 

4 

420 

59 

14.86 

1.  13 

3.  14 

5.35 

6.22 

5 

36  0 

45 

14.43 

1.56 

2.58 

4.01 

5.06 

5 

30  0 

62 

14.32 

1.70 

3.64 

5.27 

6.49 

7 

24  0 

51 

1  3.41 

2.52 

3.  36 

4.22 

5.52 

8 

180 

61 

13.27 

1.94 

3.  33 

4.49 

5.66 

9 

120 

51 

12.52 

2.48 

3.04 

4.01 

5.30 

10 

60 

53 

1  2.33 

2.67 

3.66 

4.74 

5.49 

1  1 

0 

54 

12.09 

1.99 

2.88 

4.09 

4.79 

Not  es  * 

1.  ~Buoy  column   indicates  station   number    in    vertical 

temperature   cross -sections. 

2.  Position  column   lists   the  range    (in   nautical    miles) 
from  the   most   seaward   buoy. 

3.  dT/5m    is    the  temperature   change    over    the    first    five 
meter   interval    in    the  thermccline;       dT/15   is    the 
temperature   change  over    the    first    fifteen   mater 
interval    in    the   thermae  line;      dT/50    and    dT/100 
follow    in   the    same   manner. 


TABLE    V 
Models    of    MLD   for    15    November. 

Model   Variables  Coefficient  of   Determination 

SST  0.111 

SST,    Position  0.818 

SST,    Position, 

Gradient    (100)  0.880 

Notes  : 

1.  SST  is  sea  surface  temperature. 

2.  Position  refers  to  station  position 
relative  to  the  most  sa award  station. 

3.  Gradient  (100)  refers  to  the  gradient 
computed  from  the  upper  100  meters  of 
the  thermocline. 

4.  The  coefficient  of  determination  will  equal 
one  for  a  perfect  fit  of  the  curve  to  data. 
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temperature  are  easily  measured  by  satellite.  Gradient  -hen 
becomes   the  primary   factor    of   uncertainity . 

2 •      Flight   Track    Two 

The  vertical  temperature  cross-section  for  17 
November  is  shown  in  Figure  4.6.  Twenty-three  stations  were 
used  to  construct  this  transect.  The  profile  is  interesting 
for   three  reasons: 

1.  This  flight  has  no  observational  "holes"  along  the 
track; 

2.  Satellite  imagery  from  NOAA-6  is  available  for  a 
portion   of   the   track; 

3.  The  vertical  temperature  structure  of  a 
satellite-observed  sea-surface- temperature  anomaly  can 
be   studied. 

At  the  top  of  the  thermocline,  the  gradient  appears 
strong  between  stations  13  and  16  (isotherms  are  densely 
packed).  Beginning  at  station  nine,  the  isotherms  begin  to 
diverge.  This  indicates  a  weakening  of  the  gradient  in  the 
upper  five  to  fifteen  meters  of  the  therraocline.  As  with 
flight  one,  the  initial  gradient  (from  five  meters)  is 
inversely  correlated  with  position  along  this  track.  The 
initial  strength  of  the  gradient  at  the  top  of  the  thermo- 
cline    decreases   towards   the    coast. 

The  location  of  the  divergence  of  the  isotherms  a4- 
the    top   of      the   thermocline     (station    9)  corresponds    to   the 

location  of  the  change  from  the  Central  Subarctic  Domain  to 
the   Transitional  Domain   as    presented    by   Dodimead      e_^_   al. 

Figure  4.7  is  a  plot  of  sea  surface  temperature  and 
mixed-layer-depth  variation  with  position.  From  station 
nine  to  station  twenty-one,  the  sea  surface  temperature 
increases  almost  linearly  with  position,  while  mixed  layer 
depth    fluctuates      around    sixty      meters.         Satellite      imagerv 
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Figure   a. 6 


Vertical   Temperature   Cross-section 
(17    November^  . 


57 


identified  a  warm  meander  between  station  four  and  station 
twenty-xwo  (Figures  4.8  and  4.9).  This  shows  up  as  a  rapid 
increase  in  sea  surface  temperature  over  a  (longer)  distance 
from      station      four      to   station      three.  The      sea      surface 

temperature  reaches  a  maximum  in  the  vicinity  of  station 
three,  and  then  falls  rapidly  across  the  southeast  wall  of 
■■-he  meander.  The  mixed  layer  depth  also  reaches  a  maximum 
at  station  three,  and  then  falls  rapidly  across  the 
southeast      wall.  This      correlation      between      sea      surface 

temperature  patterns  and  vertical  structure  in  the  vicinity 
of  frontal  positions  was  also  noted  by  Legeckis  and  Gordon 
(Ref.    33),    cited   earlier. 

Examinations  of  this  track  based  on  domains  and 
satellite  imagery  requires  that  the  track  be  split  into 
three  separate  sections.  The  first  section,  in  the  Central 
Subarctic      Domain,  includes      stations        thirteen      through 

sixteen.  The  Transitional  Domain  is  then  split  using  the 
satellite  imagery.  One  section  extends  from  station  nine  to 
station  twenty-one.  The  warm  meander  then  constitutes  the 
final   section   and   includes    stations   four    to   twenty-four. 

Table  VI  provides  position,  mixed  layer  depth,  sea 
surface  temperature,  and  gradients  in  the  thermociine  for 
each  station  occupied  for  flight  track  two.  As  with  flight 
one,  this  table  constitutes  the  matrix  used  to  examine 
relationships  between  various  elements  which  may  lead  to  a 
prediction    of   mixed   layer  depth. 

Correlations  of  sea  surface  temperature  and  mixed 
layer  depth  along  the  entire  track  yield  low  values  for  the 
correlation  coefficient  (r)  .  For  the  entire  track,  there  is 
little  correlation  between,  sea  surface  temperature  and  mixed 
layer  depth  (r  =  0.17).  Mixed  layer  depth  as  a  function  of 
position,  sea  surface  temperature,  and  gradient  from  the  100 
meter    thermociine      is    better    (r      =   0.34)  .  However,      along 
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Figure  4.7 


Mixed   Layer   Depth   and   Ssa    Surface   Temperature 
(17    November)  . 
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Figure   U.8        NOAA-6   Image    (17   November)     [Ref.    60] 
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Figure  4.9 


Flight   Track.    Buoy    Numbsrs    and   Meander   Location 
in   Figure  4.8. 
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TABL2    VI 
MLD,    SST   and   Gradients  ax   the  Thermocline,    17    November    1980- 


uoy 

Position 

MLD 

SST 

dT/5m 

dT/15m 

dT/50m 

dT/100m 

2 

ffl 

45 

13.76 

1.33 

2.  35 

3.40 

3.97 

3 

70 

15.80 

1.51 

2.82 

4.77 

6.53 

4 

510 

67 

14.83 

1.48 

3.  71 

5.01 

6.04 

5 

45  0 

54 

14.46 

1.51 

3.  47 

4.95 

6.11 

6 

39  0 

61 

13.97 

1.64 

2.  93 

4.09 

5.33 

7 

330 

54 

13.06 

1.69 

3.  25 

4.33 

5.70 

8 

270 

70 

12.87 

2.65 

3.  35 

4.64 

5.30 

9 

210 

64 

12.  14 

1.  19 

2.96 

4.36 

4.95 

10 

150 

58 

1  1.95 

1.67 

3.70 

4.76 

5.52 

1  1 

90 

58 

1  1.63 

2.58 

3.  39 

4.73 

5.  19 

12 

30 

58 

10.61 

1.85 

2.  63 

3.87 

4.23 

13 

0 

56 

10.13 

2.47 

3.09 

3.90 

4.25 

14 

60 

48 

10.72 

1.71 

2.68 

3.62 

4.24 

15 

120 

56 

1  1.71 

2.58 

3.  36 

4.47 

5.01 

16 

180 

58 

1  1.83 

2.03 

3.72 

4.82 

5.25 

17 

24  0 

53 

12.65 

1.75 

2.  3  8 

3.90 

4.95 

18 

300 

59 

12.90 

1.43 

2.  83 

4.07 

5.64 

19 

36  0 

53 

13.86 

2.66 

3.63 

5.06 

6.22 

20 

42  0 

54 

1  4.43 

1.  16 

3.55 

5.20 

6.24 

21 

480 

58 

14.51 

1.83 

3.  84 

4.  34 

6.01 

22 

54  0 

53 

15.24 

1.80 

3.  50 

5.70 

6.70 

23 

60  0 

69 

15.59 

0.93 

2.  64 

4.66 

6.  30 

24 

66  0 

24 

13.54 

1.08 

1.95 

4.28 

5.35 

shorter  track  segments,  and  particularly  across  the  meander, 
mixed  layer  depth  and  sea  surface  temperature  show  more 
correlation  than  for  the  entire  track.  Across  the  meander, 
the  correlation  coefficient  for  mixed  layer  depth  vs.  sea 
surface    temperature      is    0.8  6.  This    was        not    the      case    in 

either  the  Central  Subarctic  Domain  (r  =  0.40)  or  in  the 
other   part    of   the   Transitional  Domain    (r   =   -0.53). 

The  mixed  layer  depth  and  sea  surface  terapera+ure 
did  not  correlate  to  any  siginificant  degree  on  17  November. 
Therefore,  only  two  models  are  proposed  for  this  flight 
track.  As   stated      previously,         breaking     the    track      into 

segments  yields  a  much  higher  correlation  between  otherwise 
uncorrelated  parameters.  Table  711  presents  the  models  for 
the    entire    track   only. 
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TABLE    VII 
Models    of   MLD   for    17    November. 

Model   Variables  Coefficient   of   Determination 

SST,    Position  0.615 


SST,    Position, 
Gradient    (15) 


0.627 


SST.    Position. 

Gradient    (100)  0.707 

Note:  , .  ,    , 

1.      Two    gradient    models    were   computed 
because   the    gradient    from    15    meters 
appeared   to    be    more    highly   correlated 
wiuh  track   position   than   the   gradient 
computed    from   the    upper    100    meters    of 
the" thermocline. 


The  higher  correlation  between  sea  surface  tempera- 
ture and  mixed  layer  depth  in  the  warm  meander  may  provide  a 
means  of  estimating  the  depth  of  the  mixed  layer  in  such 
anomalies   through   observations  of   sea    surface   temperature. 

During  the  interval  between  flight  one  and  flight 
two,  one  storm  had  traversed  the  area  and  another  had  begun 
to      approach   the      most   western      stations    on      the    track.  A 

slight  deepening  of  the  10-degree  through  11-degree 
isotherms  occurred.  The  average  mixed  layer  depth  changed 
from  53.7  meters  on  15  November  to  57.0  meters  on  17 
November. 

3 •      Flight   Track    Three 

Figure  4. 10  is  the  vertical  cross-section  associated 
with  flight  track  three.  Fourteen  stations  were  used  to 
construct   this   transect.  Although   there      are   some    holes    in 

the  data  along  the  track,  this  profile  presents  a  picture 
similar   to    flight    track   two. 
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The  gradient  at  the  top  of  the  thermocline  diverges 
beginning  at  station  nine,  just  as  in  the  previous  sections. 
Again,  the  gradient  computed  from  the  first  five  meters  of 
the  thermocline  was  inversely  correlated  with  position  (r  = 
-0.766)  .  The  meander  discussed  for  flight  two  still  can  be 
located   in    this    profile   between    stations    four   and   two. 

The  sea  surface  temperature  and  mixed  layer  depth 
variations  with  position  are  plotted  in  Figure  4.11.  Table 
VIII  lists  the  values  for  mixed  layer  depth,  sea  surface 
temperature,  and  gradients  alcnc  the  flight  track.  Between 
station  nine  and  station  six  (Transitional  Domain)  both  sea 
surface  temperature  and  mixed  layer  depth  vary  almost 
linearly.  Over   the      entire   track      length,      however,         the 

correlation  between  the  two  is  negligible  (r  =  0.10) .  From 
station  twenty-one  to  station  one  (across  the  meander)  the 
correlation  between  mixed  layer  i^pth  ard  s»>  surface 
temperature   is   much   better     (r  =   0.51). 

A  front  moved  into  the  flight  track  area  on  17 
November,  and  became  stationary  over  the  area  on  18 
November.  The  stationary  front  caused  poor  conditions  for 
satellite  observations  due  to  clouds.  Small  changes  in  sea 
surface  temperature  patterns  could  be  expected  from  the 
front,  but  the  basic  structure  remained  unchanged  from  17 
November.  The  average  mixed  layer  depth  for  19  November 
remained  at  57.0  meters.  The  position  and  structure  of  the 
meander  along  the  track  remained  unchanged  except  for  -he 
development  of  a  peak  temperature  exceeding  16  degrees 
Celsius   at   station    three. 

Table  IX  summarizes  the  models  examined  for  flight 
track  three.  The  correlation  coefficient  computed  for  sea 
surface  temperature  versus  mixed  layer  depth  was  extremely 
low  (r  =  0.01).  Therefore,  the  msdel  eguation  for  mixed 
layer    depth      as    a      function    of      sea    surface      temperature   was 
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Vertical  Temperature   Cross-section 
(19   November) . 
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Figure  4.11 


Mixed  Layer  Depth  and  Sea  Surface  Temperature 
(19  November)  . 
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TABLE   VIII 
MLD,    SST   and  Gradients  at   the  Thermos line,    19    November    1980 


Buoy  Position 


1 

48  0 

3 

36  0 

4 

300 

5 

24  0 

6 

180 

7 

120 

8 

60 

9 

0 

17 

30 

18 

90 

20 

210 

21 

27  0 

22 

33  0 

23 

39  0 

24 

45  0 

MLD 

SST 

dT/5m 

dT/15ra 

dT/50m 

dT/100m 

30 

13.17 

1.15 

2.86 

3.90 

4.65 

67 

16.04 

1.41 

2.69 

4.62 

5.92 

58 

14.81 

1.  13 

3.  14 

5.35 

6.22 

59 

14.43 

1.56 

2.58 

4.01 

5.06 

54 

13.95 

1.70 

3.64 

5.27 

6.49 

53 

13.14 

2.52 

3.  36 

4.22 

5.52 

62 

1  3.00 

1.94 

3.33 

4.49 

5.66 

66 

1  1.98 

2.48 

3.04 

4.01 

5.30 

64 

12.79 

2.67 

3.66 

4.74 

5.49 

61 

12.95 

1.99 

2.88 

4.09 

4.79 

64 

14.32 

1.  15 

2.  86 

3.90 

4.65 

54 

14.35 

1.41 

2.  69 

4.62 

5.92 

43 

15.80 

1.  13 

3.  14 

5.35 

6.22 

59 

15.48 

1.56 

2.58 

4.01 

5.06 

54 

13.51 

1.70 

3.  64 

5.27 

6.49 

omitted.  Instead,  a  second  gradient  model  was  introduced  as 
with  flight  track  two.  The  gradient  from  100  meters  was 
continued  as  an  element  in  a  separate  model  because  this 
gradient  was  either  the  most  influential  or  second  most 
influential   of   all    gradients   used   in    model    equations. 

The  best  fit  of  any  model  for  mixed  layer  depth  was 
the  model  based  only  en  sea  surface  temperature  and  track 
position.  Both  surface  temperature  and  track  position  ar? 
variables  which  can  be  derived  from  a  satellite  image  of  the 
sea    surface. 

*«      Flight   Track    Four 

The  vertical  cross- section  of  temperature  for  flight 
track  four  (Figure  4.12)  was  constructed  from  21  stations. 
It  is  the  most  complete  transect  except  for  17  November.  A 
definite  cooling  has  occurred  between  17  November  and  this  1 
December  -rack.  Surface  temperaturas  dropped  by  about  one 
degree  Celsius  (-.he  average  sea  surface  temperature  dropped 
from    13.1      degrees    Celsius      on    17      November   to      12.0    degrees 
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TABLE    IX 
Models    of    SLD   for    19    November. 

Model   Variables  Coefficient   of  Determination 

SST,    Position  0.615 

SSI t    Position,  „   _    „ 

Gradient    (5)  0.547 

SST.    Position, 

GraSien"     (100)  0.517 

Note: 

1.      Two   gradient   models    were   compute! 
because  the    gradient    from  5    meters 
appeared   to   be    more   highly  correlated 
with  track   position  than  the   gradient 
computed   from   the    upper    100    meters   of 
the  thermocline. 


Celsius  on  1  December).  The  gradient  in  the  upper  five  to 
fifteen  meters  of  the  thermocline  is  not  as  strong  as  along 
the      November      flight    tracks.  Although      still      negatively 

correlated  with  track  position,  the  correlation  between 
gradient  and  track  position  is  much  weaker  for  this  track 
(r  =  -0.16)  than  for  similar  gradients  computed  from  the 
upper  five  meters  of  the  thermocline.  The  meander  is  still 
recognizable,  but  its  vertical  profile  is  much  less 
pronounced   than   during  the    November   flights. 

The  gradient  computed  from  the  upner  five  meters  of 
the  thermocline  showed  the  strongest  correlation  with  posi- 
tion for  all  gradients  computed  along  flight  track  four.  As 
noted  previously,  the  correlation  between  the  five  meter 
gradient  and  position  has  steadily  decreased  since  the 
first  flight.  However,  the  five  meter  gradient  is  the  most 
highly  correlated  of  any  gradient  with  position  on  each  of 
the    four  flights. 

Mixed  layer  depth  doesn't  show  much  relationship  to 
track    position.       Sea    surface   temperature    dees    increase    along 
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the    track,         but   the      change   across      the    meander      is    not  as 

significant   as   in      past    flights.        Figure    4.  13      displays  the 

variation   with   position     for    both   mixed    layer      depth    and  sea 
surface   temperature. 


TABLE    X 
MID,    SST   and   Gradients  at    the  Thermocline,    1    December    1980 


iuoy 

position 

HLD 

SST 

dT/5m 

dT/15m 

dT/50m 

dT/100m 

1 

69  0 

51 

12.98 

1.  15 

2.86 

3.90 

4.65 

4 

510 

61 

14.46 

1.41 

2.69 

4.62 

5.92 

5 

45  0 

64 

13.30 

1.13 

3.  14 

5.35 

6.22 

6 

39  0 

69 

12.92 

1.56 

2.58 

4.01 

5.06 

7 

330 

72 

1  1.74 

1.70 

3.64 

5.27 

6.49 

8 

27  0 

72 

1  1.82 

1.94 

3.33 

4.49 

5.66 

9 

210 

61 

1  1.  12 

2.52 

3.  36 

4.22 

5.52 

10 

150 

72 

10.72 

1.94 

3.  33 

4.49 

5.66 

1  1 

90 

78 

10.77 

2.52 

3.  36 

4.22 

5.52 

12 

30 

70 

9.13 

2.48 

3.04 

4.01 

5.30 

13 

0 

74 

8.94 

2.67 

3.66 

4.74 

5.49 

14 

60 

72 

9.51 

1.99 

2.88 

4.09 

4.79 

15 

120 

69 

10.32 

1.  15 

2.36 

3.90 

4.65 

17 

24  0 

70 

1  1.42 

1.41 

2.69 

4.62 

5.92 

18 

300 

78 

1  1.70 

1.  13 

3.  14 

5.35 

6.22 

19 

36  0 

74 

1  2.98 

1.56 

2.58 

4.01 

5.06 

20 

420 

80 

13.25 

1.70 

3.64 

5.27 

6.49 

21 

48  0 

69 

13.54 

1.99 

2.  38 

4.09 

4.79 

22 

540 

66 

1  4.24 

2.6^ 

3.66 

4.74 

5.49 

23 

600 

59 

14.27 

2.48 

3.04 

4.01 

5.30 

24 

66  0 

50 

13.06 

2.67 

3.  66 

4.74 

5.49 

Table  X  lists  the  factors  which  were  examined  for 
flight  four.  This  flight  was  unusual  in  that  large  changes 
from    previous   flights    occurred    in    tha    correlations. 

Both  sea  surface  temperature  and  mixed  layer  depth 
had  higher  correlations  with  track  position  than  on  previous 
flights.  The      correlation     of        gradients      with      position 

declined  in  all  four  cases  examined.  Mixed  layer  depth  was 
more  highly  correlated  with  sea  surface  temperature  than  on 
previous   flights. 
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Figure  4.12 


Vertical  Temperature   Cross-section 
(1    December)  . 
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Pigure  4.13 


Mixed  Layer  Depth  and  Sea  Surface  Temperature 
(1  December)  . 
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Three      models      ar  =      presented    for      flight      four      and 
summarized      in   Table      XI.  Since      the   correlation      of      sea 

surface  temperature  and  mixed  layer  depth  was  the  highest 
observed  of  all  flights,  i t  is  included  as  a  model.  Mixed 
layer  depth  as  a  function.  of  track  position,  sea  surface 
temperature  and  thermocline  gradient  showed  the  best 
correlation  when  the  gradient  computed  from  100  meters  was 
used. 


TABLE    XI 
Models    of   MLD    for    1    December. 

Model   Variables  Coefficient   of   Determination 

SST  0.175 

SST,    Position  0.6  09 

SST.    Position, 

Gradient    (100)  0.663 

Note: 

1.      One   gradient   model  was    computed 

because  the    aradient    from   100   meters 
appeared   to    Be    more   highly   correlated 
with   track   position   than  any    other 
gradient. 


5«      Flight    Track    Five    (Center   and    South) 

The  fifth  flight  added  another  dimension  to  ■'"he  AXBT 
analysis.  The  southern  track,  offset  by  60  nautical  miles 
from  the  center  track,  provided  additional  information  about 
the  extent  of  the  warm  meander  discussed  previously.  Figure 
4.14,  constructed  from  eleven  stations,  is  the  vertical 
Temperature  Cross-section  for  the  center  track  on  3 
December.         Figure      4.15,      constructed    from      eight    sta-ions, 
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shows  the  profile  for  the  southern  track.  Table  XII  lists 
the  values  for  mixed  layer  depth,  sea  surface  temperature, 
and    gradients      along    the   center      flight    track.  Table   XIII 

lists   similar   data    for  the    southern   track. 

The  vertical  temperature  profiles  for  the  two  tracks 
are   similar      in   some    respects.  Ih=    meander      is    definitely 

present  en  the  center  track.  The  southern  track  displays 
some  similar  subsurface  structure  in  the  area  where  the 
meander  should  be  located  (between  stations  forty-one  and 
forty-three).  Although  station  forty-two  should  be  inside 
the  warm  meander,  the  surface  temperature  at  the  station  is 
13.5  degrees  Celsius,  about  one-half  degree  too  low  to  be  a 
part  of  the  meander.  The  closed  warm  pocket  of  water  is  not 
present  in  the  southern  track  (the  center  track  has  a  closed 
pocket   of    14    degree   Celsius    water). 

The  meander/eddy  influence  appears  to  be  present  to 
about  290  meters  in  both  tracks.  The  southern  track  does 
not  have  the  pronounced  structure  shown  along  the  center 
track.  Although  some  signs  of  ths  meander  remain  in  the 
southern  track,  the  meander  appears  to  have  formed  a  closed 
circulation. 

Along  the  center  flight  track,  sea  surface 
temperature  and  mixed  layer  depth  (Figure  4.16)  are  not 
highly  correlated  (r  =  0.38).  Across  the  eddy,  however,  the 
correlation  coefficient  for  sea  surface  temperature  versus 
mixed  layer  dep+h  increases  to  0.88.  The  southern  track 
(Figure  4. 17)  also  shows  an  excellent  correlation  between 
sea  surface  temperature  and  mixed  l=yer  depth  (  r  =  0.84). 
The  corrleation  between  these  two  variables  is  less  (r  = 
0.66)  across  the  track  segment  which  coincides  with  the 
eddy.  It  appears  that  the  meander  has  formed  into  an  eddy 
which    dees   not    extend    to    the   southern    track. 
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Figure   4.14 


Vertical  Temperature   Cross-section 
(3   Dec,    Center)  . 
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Figure   U.15 


Vertical  Temperature   Cross-section 
(3    Dec,    South) . 
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Figure  a.  16 


Mixed  Layer  Depth  and  Sea  Surface  Temperature 
(3  Dec,  Center)  . 
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Figure  4.17 


Mixed  Layer  Depth  and  Sea  Surface  Temperature 
(3  Dec,  South)  . 
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Most  of  the  southern  track  is  in  the  Transitional 
Domain.  The  Central  Subarctic  Domain  has  influence  up  to 
station  thirty-five,  after  which  the  track  is  under  the 
influence  of  the  Transitional  Domain.  Changes  apppear  to  be 
linear   along   the   track  within  the   Transitional   Domain. 

Winds  increased  from  11  meters/sec  (24  kts)  to  14-16 
meters/sec  (30-35  kts)  on  1  December.  These  wind  speeds 
were  sustained  through  2  December  and  through  most  of  3 
December  [Ref.  61].  Thess  cold  northwesterly  winds  lowered 
the  sea  surface  temperature  along  the  track,  but  caused  no 
major   changes  to  the    vertical  Temperature   Cross-section. 

The  average  mixed  layer  depth  for  the  center  track 
was  66.4  meters  on  3  December.  The  southern  track  averaged 
64.6    meters.  There    was      a   decrease      in   the      average   mixed 

layer  depth  along  the  track  since  the  previous  flight  on  1 
December.  In  comparing  the  two  days,  a  large  decrease  in 
mixed  layer  occurred  at  station  five  (21  meter  change).  At 
first  glance  this  appears  to  be  an  instrument  error. 
However,  the  mixed  layer  depth  profile  for  flight  track  one 
(Figure  4.5)  is  very  similar  to  the  profile  for  flight  track 
five.  On  both  days  station  five  had  a  shallow  mixed  layer 
which    deepened   considerably    across   the    eddy. 

When  overlaid,  the  mixed  layer  depth  profiles  and 
the  sea  surface  temperature  profiles  for  tracks  one  and  five 
are    remarkably   similar.  The   mixed    layer   depth      for   flight 

five  is  deeper  than  for  flight  one,  but  the  variations  with 
position  are  very  similar.  Between  stations  ten  and  two, 
the  sea  surface  temperature  is  highly  correlated  (r  =  0.98) 
with  pcsition,  and  the  mixed  layer  dep^h  shows  some 
correlation    (r   =   0.54). 

Four  model  equations  were  used  for  each  track  on 
flight  five.  The  models  are  summarized  in  Table  XIV  for  the 
center  track,      and    in   Table    XV      for   the    southern    track.        On 
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TABLE  XII 
MLD ,  SST  and  Gradients  at  the  Thermccline,  3  December  1980, 


Buoy  Position 

46(1)       69  0 

3  570 

4  510 

5  45  0 

6  390 

7  33  0 

8  270 

9  210 
10  150 

12  30 

13  0 


MLD 

ssr 

dT/5m 

dT/15m 

dT/50m 

dT/100m 

45 

12.47 

0.04 

1.  28 

3.58 

4.48 

70 

14.64 

1.05 

2.  39 

4.67 

5.97 

62 

14.16 

1.35 

2.73 

4.87 

5.46 

42 

13.38 

1.41 

2.  17 

3.35 

4.51 

72 

13.06 

1.42 

2.44 

4.20 

5.20 

69 

1  1.77 

0.81 

2.07 

3.01 

4.  14 

83 

1  1.42 

1.53 

2.66 

3.71 

4.36 

77 

10.99 

1.80 

2.  26 

3.24 

3.37 

62 

10.02 

2.28 

2.71 

3.30 

3.82 

66 

8.86 

0.78 

1.  43 

2.24 

2.43 

82 

8.46 

1.42 

1.88 

2.20 

2.49 

TABLE    XIII 

MLD,  SST  and  Gradients  at  the  Thermocline,  3  December  1980 

(South) 


Buoy  Position 


34 

0 

35 

60 

33 

24  0 

39 

30  0 

41 

420 

42 

48  0 

43 

54  0 

44 

600 

MLD 

SST 

dT/5m 

dT/15m 

dT/50m 

dT/100m 

80 

9.40 

1.23 

2.  56 

3.23 

3.42 

86 

9.53 

1.43 

1.70 

2.40 

2.62 

70 

1  1.44 

1.69 

1.  38 

3.08 

3.55 

64 

12.04 

1.00 

2.  53 

3.34 

4.61 

74 

13.27 

0.92 

2.60 

3.43 

4.89 

43 

13.46 

1.79 

3.  61 

4.20 

4.76 

46 

13.89 

1.62 

2.  56 

4.20 

5.  12 

54 

13.65 

1.07 

2.  10 

4.09 

5.22 

the  south  track,  the  gradient  computed  from  100  meters 
produced  the  worst  fit  of  any  model  of  mixed  layer  as  a 
function  of  sea  surface  temperature,  track  position  and 
thermocline  gradient.  Th=  gradient  from  100  meters  produced 
the  best  fit  for  the  center  track  aodel.  For  comparison, 
the  southern  track  model  using  the  gradient  computed  from 
100  meters  is  included. 
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TABLE    XIV 


Models   of    MLD    for    3    December,    Center    Track 


Model   Variables 

SST 

SST,    Position 

SST,    Position, 
Gradient    (15) 

SST.    Position. 
Gradient    (100) 


Coefficient   of   Determination 
0.146 
0.404 

0.441 

0.523 


TABL2    XV 


Models    of    MLD  for    3   December,    South    Track 


Model   Variables 

SST 

SST,    Position 

SST,    Position, 
Gradient    (50) 

SST.    Position. 
Gradient    (100) 


Coefficient   of  Det ermination 
0.703 
0.716 

0.861 

0.716 


6-   fliakt  Track  Six  (Center  and  North) 

Flight  six  examined  both  the  center  track  (Figure 
4.18)  and  a  parallel  track  60  nautical  miles  to  the  north 
(Figure  u..  19)  with  eleven  stations  each.  The  values  of 
mixed  layer  depth,  sea  surface  temperature,  and  gradient 
associated  with  each  station  are  given  in  Table  XVI  for  the 
center  track,  and  Table  XV II  for  the  northern  track.  The 
temperature   structure  along   the   center  track   res  ambles 
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previous  flight  tracks.  The  isotherms  begin  to  diverge 
after  station  nine.  The  gradient  computed  from  the  upper 
five  meters  of  the  thermocline  is  negatively  correlated  with 
position  (r  =  -0.27). 

The  eddy  seems  to  have  diminished  both  in  length 
along  the  center  track  and  in  depth.  The  northern  track 
does  not  appear  to  traverse  the  eddy. 

The  sea  surface  temperature  and  mixed  layer  depth 
variation  with  position  along  the  center  track  (Figure  4.20) 
show  the  best  correlation  in  any  flight  along  the  center 
track.  Both  sea  surface  temperature  and  mixed  layer  depth 
varied  linearly  with  position  up  to  station  three.  Sea 
surface  temperature  correlation  with  position  along  the 
entire  track  was  better  (r  =  0.92)  than  along  the  shorter 
track  across  the  eddy  (r  =  -0.80) ,  cr  along  any  previous 
center  flight  track.  Mixed  layer  depth  was  more  highly 
correlated  with  position  across  the  eddy  (r  =  -0.95)  than 
along  the  entire  track  (r  =  -0.83). 

Across  the  eddy  on  the  center  track,  sea  surface 
temperature  and  mixed  layer  depth  again  are  highly 
correlated  (r  -  0.90)  .  Dver  the  entire  track  length,  the 
relationship  between  sea  surface  temperature  and  mixed  layer 
depth  is  not  as  strong  (r  =  -0.55).  Of  ail  flights,  this 
was  the  best  overall  correlation  of  sea  surface  temperature 
and  mixed  layer  depth  for  the  center  track. 

The  northern  track  (Figure  4.21)  shows  a 
relationship  similar  to  that  along  the  center  track  between 
sea  surface  temperature  and  mixed  layer  depth.  The  eddy  is 
not  apparent  in  the  sea  surface  temperature  profile,  but 
does  present  a  familiar  profile  of  mixed  layer  depth. 
Across  the  eddy  position  on  the  north  track,  mixed  layer 
depth  is  highly  correlated  with  position  (r  =  -0.97).  This 
same  relationship  holds  true  for  kht  center  track  (r  = 
-0.95)  . 
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Figure   4.18 


Vertical  Temperature   Cross-section 
(5   Dec,    Center)  . 
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Figure  4.19 


?ertical  Temperature   Cross-section 
(5   Dec,    NortS) . 
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Figure   4.20 


Bixed   Layer   Depth    and   Sea    Surface  Temperature 
(5   Dec,    Center)  . 
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Figure  U.21 


Mixed  Layer  Depth  and  Sea  Surface  Temperature 
(5  Dec,  North) . 
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A  subsurface  featura  developed  over  the  course  of 
tha  experiment,  and  became  very  obvious  in  the  northern 
track  on  5  December.  At  station  58,  the  isotherms  reach  a 
maximum  depth.  A  corresponding  excursion  to  depth  for  the 
same  isotherms  (6  and  7  degree  Celsius)  occurs  along  •'-he 
center  track  at   station  nine.  In  looking   back   through   the 

previous  flight  cross-sections,  this  same  excursion  is  found 
to   occur  on    flights    four   and    five. 

Lundell  [ Ref .  62]  reported  a  submerged  sound  channel 
in  the  vicinity  of  station  nine.  The  Kilonski  digitized 
bathythermographs  indicated  a  slight  inversion  in  this  area 
(of  the  order  0.3  degrees  Celsius).  This  unusual  structure 
occurs   at      the   boundary    between      two    Domains:  the   Central 

Subarctic  and  Transitional.  Along  the  center  track,  the 
maximum  depth  excursion  of  the  six  and  seven  degree 
isotherms  shifted  to  the  west  between  flight  two  and  flight 
six. 

This  submerged  feature  is  not  driven  by  surface 
conditions  as  it  occurs  well  below  both  thermocline  and 
halocline.  This  appears  to  be  some  type  of  shift  in  lower 
domain  boundaries.  A  time  series  over  one  or  two  months, 
between  stations  seven  and  ten,  might  have  provided  some 
insight  into  the  lower  layer  variability  at  this  domain 
boundary. 

As  with  the  southern  track  on  flight  five,  the 
northern  track  data  on  flight  six  showed  the  worst  fit  by  a 
model    using      the   gradient      computed    from      100    meters.  The 

gradient  from  fifty  meters  produced  the  best  fit  on  this 
track.  The  gradient  from  five  meters  produced  the  best  fit 
on  the  center  track.  The  model  parameters  and  coefficients 
of   determination   are    presented   in   Tables   XVIII   and   XIX. 
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TABLE   XVI 
MLD,    SST   and   Gradients  at   the  Thermocline,    5    December    1980 


uoy 

Pcsiti 

on 

MLD 

SST 

dT/5m 

dT/15m 

dT/50m 

dT/100m 

1 

690 
63  0 

42 

12.63 

0.84 

1.  17 

3.  19 

4.21 

2 

58 

13.11 

1.77 

3.  17 

3.79 

4.40 

3 

570 

70 

14.40 

0.89 

2.  10 

4.25 

5.88 

4 

510 

72 

13.81 

1.79 

3.05 

4.60 

5.44 

6 

39  0 

78 

12.63 

1.77 

3.  05 

4.00 

4.83 

7 

330 

75 

1  1.77 

1.  10 

2.  18 

3.40 

4.30 

8 

27  0 

78 

1  1.23 

1.39 

2.  45 

3.67 

4.25 

9 

210 

75 

1  0.88 

1.  11 

2.46 

3.42 

3.42 

10 

150 

78 

9.94 

1.76 

2.57 

3.47 

3.75 

12 

30 

82 

8.97 

1.60 

2.  26 

2.91 

2.83 

13 

0 

83 

8.43 

1.46 

1.83 

2.02 

2.32 

te: 
A 

s: 

Light 

invers 

icn   exi 

sted  at 

station 

12    bet 

ween 

135 

mete 

rs 

and 

20  0   met 

s  rs. 

TABLE  XVII 

MLD,  SST  and  Gradients  at  the  Thermocline,  5  December  1980 

(North) 


Buoy    Position      MLD       SST         dT/5m      dT/15m      dT/50m      dT/l00m 


54 
55 
56 
57 
58 
59 
60 
61 
62 
64 
65 


0 
60 
120 
18  0 
240 
300 
36  0 
42  0 
48  0 
600 
66  0 


70 
83 
74 
77 
64 
59 
70 
69 
69 
50 
45 


8.81 
9.94 
10.15 
10.37 
10.96 
1  1.20 
1  1.93 
12.60 
12.73 
12.84 
12.63 


1.22 
1.36 
2.03 
2.24 
0.91 
1.70 
0.95 
0.92 
1.64 
1.86 
1.  12 


1.59 
2.  32 


2. 
2. 
2. 

Am    • 

1. 

2. 
3.  42 
3.01 
1.  44 


68 
89 
16 
31 

97 
28 


2.  14 
3.39 
3.65 
3.79 
2.39 
3.28 
3.03 
3.43 
3.96 
3.87 
2.  28 


37 
69 
77 
23 
69 
49 
52 
89 
62 
71 
59 


E.       DYNAMIC    HEIGHT    COMPARISONS,     17    NOVEMBER 

Emery  and  O'Brien  [Ref.  63]  proposed  a  method  of  infer- 
ring salinity  from  temperature  or  depth  in  the  North  Pacific 
Ocean.        The    method    relies    on      mean   salinity   values    computed 
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TABLE    XVIII 


Models   of   MLD    for   5    December,    Center  Track 


Model   Variables 

SSI 

SST,    Position 

SST.    Position, 
Gradient    (5) 

SST,    Position, 
Gradient    (100) 


Coefficient    of  Determination 
0.298 
0.962 

0.977 
0.963 


TABLE    XIX 


Models    of    MLD   for    5    December,    North   Track. 


Model   Variables 

SST 

SST,    Position 

SST,    Position, 
Gradient    (50) 

SST,    Position. 
Gradient    (100) 


Coefficient   of   Determination 
0.370 
0.8U6 

0.372 

0.85U 


by      five   degree      squares.         In      the    test      region,       the      mean 
salinity      curves   are      nearly      constant.  Bernstein    ej:      al. 

[Ref.  6U  ]  used  this  method  to  estimate  dynamic  heights  in 
the  Western  North  Pacific  with  good  success.  It  was  decided 
to  employ  both  a  single  recent  salinity  sounding  in  the  area 
(from  a  CTD  cast) ,  and  a  Flee*  Numerical  Oceanographic 
Center  salinity  climatology  for  November,  as  a  basis  for 
calculating  and  comparing  dynamic  heights  for  flight  track 
two     (17    November)  . 

88 


TABLE  XX 
Dynamic  Heights  Relative  to  350  meters  (17  Nov) 


ODY 

POSITION 

CONST-S 

FNOC-S 

13 

38: 

0.618 

0.629 

12 

0.643 

0.6*0 

14 

60. 

0.638 

0.644 

11 

90. 

0.654 

0.658 

15 

120. 

0.654 

0.658 

10 

150. 

0.656 

0.657 

16 

180. 

0.665 

0.663 

9 

210. 

0.695 

0.690 

17 

240. 

0.705 

0.698 

8 

270. 

0.713 

0.698 

18 

300. 

0.717 

0.690 

7 

330. 

0.714 

0.681 

19 

360. 

0.7U2 

0.699 

6 

390. 

0.741 

0.761 

20 

420. 

0.743 

0.753 

5 

450. 

0.759 

0.768 

21 

480. 

0.766 

0.774 

4 

510. 

0.777 

0.779 

22 

540. 

0.763 

0.747 

3 

570. 

0.821 

0.821 

23 

600. 

0.817 

0.816 

2 

630. 

0.7  37 

0.735 

24 

660. 

0.704 

0.697 

Not  es  z 

1.  "Heights   in  dynamic  meters. 

2.  Const-S   refers   to    sinale   salinity    profile   for    entire   area 

3.  FNOC-S    refers   to    November  salini- y* climatology   from   ^.he 
Fleet    Numerical   Oceanographic  Center. 


Consideration  then  was  given  to  plotting  the  dynamic 
topography  along  the  track.  Since  only  one  salinity  profile 
was  used  to  calculate  the  dynamio  height  at  all  stations, 
the  resulting  variations  in  dynamic  heights  along  the  track 
were  due  to  temperature  variations  only.  The  reference 
level  was  selected  as  350  meters  because  this  depth  repre- 
sented the  useful  limit  of  temperature  da-a  available  from 
the    digitized    report. 

ks  a  means  of  comparing  dynamic  height  variation  with 
position  and  other  factors  of  interest  such  as  temperature 
gradient  at  the  thermocline  or  mixed  layer  depth,  the  Fleet 
Numerical   Oceanographic      Center   climatology      for    the      flight 
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track,   was      used   in      calculations   of      dynamic   heights.  The 

results   are   shown   in    Figure    u.22    and   Table   XX. 

The  variation  of  dynamic  height  along  the  track  was  low 
(0.19  meters/kilometer).  Across  the  eddy  the  height  change 
was  0.1  meters  per  150  kilometers.  This  height  signal  would 
be  too  small  for  the  SEASAT  aitimetsr  to  detect  accurately. 
However,  future  generations  of  satellite  altimeters  should 
be  able  to  detect  the  magnitude  of  the  signal  from  this 
eddy. 

The  temperature  gradient  in  the  ^hermociine  can  be  an 
important  factor  in  a  model  of  mixed  layer  depth.  This 
gradient  cannot  be  predicted  accurately  from  sea  surface 
temperature  observations.  Therefora,  a  correlation  between 
thermocline  gradient  and  dynamic  height  was  investigated. 
The  gradient  from  one  hundred  meters  showed  the  strongest 
correlation  with  dynamic  height.  The  possibility  exists, 
then,  to  predict  the  gradient  from  measured  values  of 
dynamic    heights    along    the   track. 

F.       OVERALL    OBSERVATIONS 

1 .      An    Ed^l   on    the  S out heas tern    Part    of  the   Track 

A    warm      eddy    formed      in    the      observation    area.  It 

began  as  an  intrusion  of  warm  water  (the  meander  imaged  by 
satellite  on  17  November)  and  eventually  appeared  to  form  a 
closed  circulation.  Based  on  the  November  image,  the  eddy 
appeared  to  have  a  90-naut ical-mile  (165  kilometer)  major 
axis  and  a  60-nautical-mile  (110  kilometer)  minor  axis.  By 
the  end  of  the  experiment,  the  major  axis  had  decreased  to 
about  60  nautical  miles,  while  the  minor  axis  did  not 
decrease   significantly. 

The  sea  surface  temperature  and  mixed  layer  depth 
were      highly   correlated      across      th=      eddy.         The      southeast 
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Surface  Heights  (dynamic  meters) 


Figure  4.22 


Dynamic  Heights  Along 
(17  Nov)  0/350. 


he  Center  Track 
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"wall"  of  this  eddy  displayed  sharp  gradients  of  both  sea 
surface  temperature  and  mixed  layer  depth.  The  correlation 
was  consistently  strong  enough  to  infer  that  the  influence 
of  the  eddy  was  confined  to  the  vertical  area  directly  below 
the  surface  thermal  manifistation,  that  is,  the  front  which 
existed  between  the  eddy  water  and  the  surrounding  water  had 
a  slope  which  approached  ninety  degrees,  especially  to  the 
southeast  of  the  eddy. 

2-   Gradients  in  the  Thermocline 

The  gradient  computed  from  the  first  five  meters  of 
the  thermocline  showed  a  negative  correlation  with  position 
in  seven  of  eight  flight  tracks.  Th=  gradient  was  strongest 
to  the  east  of  station  nine.  Station  nine  lies  near  the 
boundary  between  the  Central  Subarctic  Domain  and  the 
Transitional  Domain. 

The  gradient  computed  from  five  meters  was  more 
correlated  with  position  in  November  than  in  December. 
Table  XXI  lists  the  correlation  coefficients  for  the 
gradient  from  five  meters. 

Of  the  gradients  computed,  the  gradient  from  five 
meters  and  the  gradient  from  one  hundred  meters  generally 
had  the  best  correlation  with  position.  Attempts  to  model 
mixed  layer  depth  were  generally  more  successful  using  the 
gradient  computed  from  100  meters.  This  is  to  be  expected 
as  the  gradient  computed  from  100  meters  tends  to  smooth 
errors  introduced  by  an  arbitrary  starting  poin+  (the 
strength  of  the  five  meter  gradient  could  vary  by  as  much  as 
twenty  percent  depending  on  what  depth  was  selected  as  the 
top  of  the  thermocline.  Table  XXII  lists  the  two  gradients, 
by  flight  track,  which  exhibited  -he  largest  magnitude 
correlation  with  position  on  the  flight  track. 
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TABLE    XXI 

Correlation  Coefficient     (r)     for   Gradient    from    5   meters 

vs. Position 


FLIGHT  DATE  CORRELATION    COEFFICIENT     (r) 

1  15    Nov  -0.818 

2  17    Nov  -0.583 

3  19    Nov  -0.766 

4  1    Dec  -0. 159 

5  3    Dec  -0.47  8 

5  (south)  3    Dec  +0.012 

6  5    Dec  -0.  266 
6  (north)            5    Dec  -0.  129 


TABLE   XXII 
Gradients   Most   Correlated  with   Position    (by   Flight) 


FLIGHT 

DATE 

GRADIENT 

<r) 

1 

15 

Nov 

5 
15 

-0.318 

-0.4  12 

2 

17 

Nov 

5 
100 

-0.583 
+0.575 

3 

19 

Nov 

5 
50 

-0.766 
+0.282 

4 

1 

Dec 

5 

50 

-0.159 
+0.115 

5 

3 

Dec 

100 
5 

+0.342 
-0.478 

5 

(S) 

3 

Dec 

100 
50 

+0.919 

+0.359 

6 

5 

Dec 

100 
50 

+0.781 
+0.542 

6 

(N) 

5 

Dec 

100 
5 

+0.524 
-0. 129 

Note: 

1.      Track    five   south    indicated   by   5(S)     and 
track   six   north   indicated  by   6  (N)  . 
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3 •      Se a   Sur f ace   Temperature 

The  sea  surface  temperature  correlated  with  position 
on  all  flight  tracks.  This  is  due  in  part  to  the  system  of 
specifying    positions   on   the    line.  It    is   still    interesting 

to  note  that  the  sea  surface  temperature  increased  almost 
linearly  along  the  track  until  reaching  the  eddy.  There,  a 
second  order  fit  with  position  was  more  appropriate.  Table 
XXIII  lists  the  correlation  coefficient  computed  for  sea 
surface   temperature   vs.    position   for   each   track. 


TABLE   XXIII 
Correlation   Coefficient    (r)    for   SST    vs.    Track   Position, 

FLIGHT  DATE  ENTIRE    TRACK         EDDY    LOCATION 

1  15    Nov  0.356  -0. 125 

2  17    Nov  0.909  -0.416 

3  19    Nov  0.569  -0.480 

4  1    Dec  0.909  -0.648 

5  3    Dec  0.894  -0.476 

5  (S)  3    Dec  0.982  0.764 

6  5    Dec  0.9  17  -0.795 
6     (N)               5    Dec                   0.954  0.491 

Not  e : 

1.      Eddy   location    is    the   along   track    correlation 
between   buoys    that   would    define   the    eddy 
boundaries    (stations    4    to   24    on    17    November) . 


The   objective    of   this   study      was    not    the    formulation 
of   a    predictor      for   sea   surface    temperature.  Observed   sea 

surface  temperature  patterns  did  shew  when  deviations  from 
the  expected  climatology  occurred.  The  sea  surface  tempera- 
ture gradients,  in  particular,  provided  very  good  informa- 
tion for  defining  frontal  conditions  at  the  eddy.  Although 
values  of  sea  surface  temperature  did  not  imply 
corresponding  values  of  mixed  layer  depth  with  any 
reliability,         the    sharp      gradients      of    surface      temperature 
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along      the      track      did  occur     in      conjunction      with      similar 
gradients   of    mixed   layer   dspth. 

Satellite-derived  sea  surface  temperature  gradients, 
then,  can  be  used  to  modify  sound  velocity  profiles  when 
water  mass  discontinuities  are  observed.  As  a  first  guess, 
a  single  bathythermograph  could  provide  temperatures,  while 
sat  el lite- derived  surface  temperatures  could  be  used  to 
modify  the  measured  vertical  structure  at  another  location. 
The  sea  surface  temperature  gradients  would  be  useful  as  a 
decision  aid  for  predicting  the  changes  occurring  in  the 
mixed  layer  depth.  A  vertical  temperature  structure  along 
the    proposed   track    could   be    worked   out    in   this   manner. 


TABLE    XXIV 

Correlation  Coefficient     (r)    -  MLD   vs.    Track   Position. 

FLIGHT  DATE  ENTIRE     TRACK         EDDY    LOCATION 

1  15    Nov  -0. 1 49  -0.231 

2  17    Nov  -0. 164  -0.607 

3  19    Nov  -0.6U1  -0.544 

4  1    Dec  -0.655  -0.931 

5  3    Dec  -0.569  -0.018 

5  (S)  3    Dec  -0.845  -0.527 

6  5    Dec  -0.327  -0.954 
6     (N)               5    Dec                -0.787  -0.972 

Note:   Eddy  location  is  the  same  as  in  the  previous  table 


**  •   Mixed  Lay_er  Depth 

The  mixed  layer  depth  varied  along  each  track. 
Sometimes  it  would  correlate  with  position,  but  seldom  would 
the  correlation  be  as  high  as  for  sea  surface  temperature. 
The  correlation  was  always  negative,  that  is  the  mixed  layer 
shallowed  toward  the  coast.   This  is  to  be  expected  from  the 
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climatological   description      of   the   area.  The   Transitional 

Domain  is  characterized  by  a  weakening  of  the  halocliner 
warmer  water,    and  shallower    isothermal    surface  layers. 

The  sea  surface  temperature  indicated  the  presence 
of  a  warm  meander  or  eddy.  The  satellite  image  from  17 
November  shows  the  thermal  signature  of  warm  water  entrained 
in  a  cold  area.  Even  so,  the  sea  surface  temperature 
generally  did  not  show  as  strong  a  correlation  with  position 
in  the  eddy  as  did  mixed  layer  depth.  Table  XXIV  shows  the 
correlation  coefficients  for  mixed  layer  depth  vs.  track 
position,  both  for  the  entire  track  and  the  eddy  only.  This 
may  provide  a  means  of  estimating  mixed  layer  depth  in  warm 
eddies.  Satellites    can      be      used      to   define     the      surface 

boundaries  of  a  warm  eddy.  Then  mixed  layer  depth  across 
the  eddy  could  be  predicted  as  a  function  of  position  in  the 
eddy. 


TABLE   XXV 
Correlation   Coefficients    (r)     for    MLD   vs.    SST. 

Flight  r    (entire   track)  r    (eddy) 

0.970 

0.858 

0.509 

0.601 

0.881 
-0.562 

0.399 
-0.433 

Note:      The    eddy   column   refers   to   correlations    of    data 

gathered    only    from   stations   along    the   track   which 
include   the    eddy. 


1 

0.334 

2 

0.170 

3 

-0.101 

4 

-0.418 

5 

-0.382 

5 

(S) 

-0.338 

6 

-0.546 

6 

(N) 

-0.608 
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5 •      Mi xed      Lay,er      Dept h      and        Sea      Surface      Temperature 
Correlations 

Mixed  layer  depth  and  sea  surface  temperature  showed 
little  correlation  in  November,  but  the  December  flight 
tracks  showed  a  strong  increase  in  the  correlation 
coefficient  for  the  entire  track.  The  eddy  displayed  a 
strong  mixed- layer- depth  sea-surface-temperature  correlation 
in    November,        decreasing   somewhat    in    December.  Table   XXV 

summarizes  the  correlation  coefficients  for  mixed  layer 
depth  and  sea  surface  temperature  for  each  day  of 
observations. 
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V.       CONCLUSIONS     AND    RECOMMENDATIONS 

A.       CONCLDSIONS 

Satellites  can  be  used  to  map  sea  surface  conditions, 
particularly     sea    surface      temperature.  For   ASTREX      data, 

knowledge  of  sea  surface  temperature  and  position  along  the 
flight  track  was  used  to  model  mixed  layer  depth.  This 
model  correlated  poorly  with  in  situ  observations.  When 
gradient  in  the  thermocline  was  added  to  the  model,  the 
correlation   increased   significantly. 

The  gradient  in  the  thermocline  may  be  predicted  (with 
low  confidence)  from  satellite  altimeter  measurements  of  the 
ocean    surface.  A   useful    relationship      between    thermocline 

gradient   and   the    dynamic   topography   variations  may   exist. 

The  most  striking  result  of  this  study  was  the  strong 
correlation  observed  between  sea  surface  temperature  and 
mixed   layer      depth    across   a      warm   eddy.  Satellite    imagery 

provided  the  location  of  the  eddy  boundaries  with  a  high 
degree   of   accuracy.  The    strength   of   this      surface   thermal 

front  is  an  indication  of  mixed  layer  depth  variability  in 
the  area.  Qualitatively,  a  large  thermal  gradient  at  the 
surface  indicates  a  correspondingly  large  change  in  mixed 
layer  depth.  Across  the  southeastern  boundary  of  the  warm 
eddy,  the  sea  surface  temperature  and  mixed  layer  depth 
changed  siginif icant ly  in  the  same  sense  (both  increased 
together    when   tracked    from    the   northwest   to   the    southeast). 

Uncorrected  satellite  temperatures  in  the  area  displayed 
the  same  variations  along  the  track  as  did  observed  AX3T  sea 
surface  temperatures.  The  variations  cf  sea  surface  temper- 
ature (thermal  gradients  at  the  surface)  contain  the  infor- 
mation     needed   to      locate   anomalies      such    as      the    warm      eddy 
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observed  in  the  south  eastern  part  of  the  flight  tracks.  The 
NOAA-6/7  satellites  can  resolve  0.5  degree  Kelvin  tempera- 
tures at  one-kilometer  resolution.  Frontal  areas  can  be 
identified  and  accurately  located.  The  AXBT's  provided 
temperature  values  at  55-kilometer  intervals,  which  were 
sufficient  to  observe  that  the  eddy  did  exist,  but  were  not 
sufficient  to  pinpoint  the  eddy  boundaries.  Satellite 
imagery,  providing  temperature  values  at  one-kilometer 
intervals,  can  be  used  to  specify  the  boundaries  of  a  front 
or   eddy   more   accurately. 

B.       RECOMMENDATIONS 

In        future        experiments        in  the        ASTREX        region, 

consideration  should  be  given  to  investigating  satellite 
located- anomalies  such  as  the  warm  eddy  observed  during 
ASTREX.  Until   a      satellite      altimeter    becomes      available, 

dynamic  heights  must  be  calculated  indirectly  from  other 
observed  variables.  Therefore,  the  salinity  should  be 
measured  in  any    future  test    area. 

Various  models  of  mixed  layer  depth  must  be  explored. 
In  addition  to  modeling  the  mixed  layer  as  a  function  of 
position,  sea  surface  temperature  and  gradient  in  the 
thermociine,      additional      factors   should      be    added.  These 

factors  should  include  wind  mixing  in  the  mixed  layer  based 
on  observed  wind  speeds  in  the  area,  dynamic  heights 
associated  with  tracks  examined,  and  possibly  observed  sea 
state   in   the   area. 

Another  model  which  could  prove  very  useful  would 
include  the  effects  of  wind  speed  and  duration  on  the  mixed 
layer  depth.  Wind  was  not  used  as  a  factor  in  this  study, 
but  it  is  an  important  element  in  the  formation  of  the  mixed 
layer.  Between  the  November  and  December  flights,  wind 
speed    increased    and   so   did      mixed    layer    depth.         Instruments 
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such  as  the  SMMR  on  NIMBU3-7  can  provide  scaler  wind  fields 
for  an  area.  The  formulation  of  mixed  layer  depth  as  a 
function  of  magnitude  and  duration  of  wind,  sea  surface 
temperatures,  and  other  factcrs  may  provide  the  best,  model 
of   mixed   layer   depth. 
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APPENDIX    A 

VARIATION    OF   THERMOCLINE    GR ADIENTS    WITH    POSITION    FOE 

NOVEMBER 
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Figure   A.  1        Gradient   fro«    5   meters   vs.    Position    (15    Nov). 
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Figure   A.  2        Gradient   froi    15   meters    vs.    Position    (15    Nov), 
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Figure    A.  3        Gradient   froa    50   meters    vs.    Position    (15    Nov). 


103 


Gradient 

(degrees   Celsius/meter) 


0.0650+ 


0.0600 


0.0550+ 


0.  0U50 


0.  ~150.  300.  U50.  600. 

TRACK   POSITION    (Nautical   Miles    from   Origin) 


Figure   A. 4        Gradient    from    100   meters    vs.  Position      (15    Nov). 
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Figure   A. 5        Gradient   fro«   5    meters    vs.  Position    (17    Nov). 
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Figure  A.  6   Gradient  from  15  meters  vs.  Position  (17  Nov). 
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Figure    &•  7        Gradient   from    50   meters    vs.  Position    (17    Nov). 
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Figure   A.  8        Gradient   from    100    meters    vs.  Position    (17    Nov). 
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Figure  A.  9   Gradient  froi  5  meters  vs.  Position  (19  Nov). 
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Figure  A. 10   Gradient  from  15  meters  vs.  Position   (19  Nov). 
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Figure   A. 11        Gradient  from    50   meters    vs.  Position      (19    Nov), 
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Figure  A. 12       Gradient  from    100    meters    vs.  Position      (19   Nov). 
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APPENDIX   B 

VARIATION    OF   THERMOCLIN  E    GR ADIENTS    WITH    POSITION    FOR 

DECEMBER 
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Figure   B.  13        Gradient   ftom  5   meters    vs.  Position      (1    Dec). 
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Figure   B.  14        Gradient   from   15    meters    vs.  Position      (1    Dec). 
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Figure   B.  15        Gradient   fron   50    meters    ts.  Position      (1    Dec). 
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Figure  B.16        Gradient  froa    100    meters    vs.  Position      (1    Dec). 
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Figure  B.  17   Gradient  fro«  5  meters  vs.  Position   (3  Dec), 
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Figure  B.  18   Gradient  froa  15  meters  vs.  Position   (3  Dec). 
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Figure  B.  19   Gradient  froa  50  meters  vs.  Position   (3  Dec). 
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Figure  B.20   Gradient  fron  100  meters  vs. Position   (3  Dec). 
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Pigure  B.21   Gradient  froa  5  meters  7s  Position   (3  Dec 
South). 
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Figure  B.22   Gradient  from  15  meters  vs  Position   (3  Dec 
South). 
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Figure   B.23        Gradient   froi    50   meters    vs   Position      (3    Dec 
South)* 
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Figure   B.24        Gradient   froa    100    meters   tts   Position       (3    Dec 
-    South). 
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Figure   B.  25        Gradient   froa   5   meters    vs.  Position      (5   Dec). 
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Figure   B- 26        Gradient   fro«    15    inetars    vs.  Position      (5    Dec), 
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Figure  B.  27   Gradient  froa  50  meters  vs.  Position   (5  Dec), 
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Figure  B.28   Gradient  froa  100  meters  vs.  Position   (5  Dec). 
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Figure   B.29        Gradient  froa    5   meters   vs.  Position       (5    Dec   - 
North). 
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Figure  B.30   Gradient  from  15  meters  vs. Position   (5  Dec  - 
North). 
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Figure  B.31        Gradient  from    50   meters    vs. Position      (5   Dec   - 
North). 
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Figure  B.32   Gradient  from  100  meters  vs.  Position   (5  Dec  - 
North). 
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AP PENDIX    C 
TABLES    OF    NORMALIZED    GRADIENTS    IN    THE    THERMOCLINE 

TABLE    XXVI 
Normalized   Gradients    for   Flight    Track   One. 

5    meter  15    meter        50    meter        100    meter 

0.1920  0.0818  0.0479 

0.2440  0.0948  0.0549 

0.2027  0.0802  0.0530 

0.2220  0.0898  0.0566 

0.2240  0.0844  0.0552 

0.2427  0.  1054  0.0649 

0.1720  0.0802  0.0506 

0.2093  0.  1070  0.0622 

0.1793  0.0924  0.0592 

0.1907  0.0780  0.0465 

Not  e  * 

1.    Gradients    were    computed    from   -he   ^emoerature 
change   over   the   distanc9    specified    in   the 
column   heading. 
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TABLE  XXVII 


Normalized  Gradients  for  Plight  Track  Two, 


5    meter 

15   meter 

50 

meter 

100    mete 

0.4940 

0.2060 

0. 

0780 

0.0425 

0.3700 

0.1753 

0. 

0774 

0.0423 

0.3420 

0. 1787 

0. 

07  24 

0.0424 

0.5160 

0.2260 

0. 

0946 

0.0519 

0.5160 

0.2240 

0. 

0894 

0.0501 

0.3340 

0.2467 

0. 

0952 

0.0552 

0.4160 

0.2480 

0. 

0964 

0.0525 

0.2380 

0.1973 

0. 

0872 

0.0495 

0.3500 

0.  1920 

0. 

0780 

0.0495 

0.5300 

0.2233 

0. 

0928 

0.0530 

0.2860 

0.1920 

0. 

0814 

0.0564 

0.3380 

0.2167 

0. 

0866 

0.0570 

0.5320 

0.2420 

0. 

1012 

0.0622 

0.3280 

0.1953 

0. 

0818 

0.0533 

0.2320 

0.2367 

0. 

1040 

0.0624 

0.3020 

0.2313 

0. 

09  90 

0.0611 

0.3660 

0.2560 

0. 

0968 

0.0601 

0.2960 

0.2473 

0. 

1002 

0.0604 

0.3600 

0.2333 

0. 

1140 

0.0670 

0.3020 

0.1880 

0. 

0954 

0.0653 

0.1860 

0.1760 

0. 

0932 

0.0630 

0.2660 

0.  1567 

0. 

0680 

0.0397 

0.2160 

0.1300 

0. 

0856 

0.0535 

TABLE  XXVIII 


Normalized  Gradients  for  Flight  Track  Three. 


5    meter 

15   meter 

50    meter 

100    mete 

0.4960 

0.2027 

0.  0802 

0.0530 

0.5340 

0.2440 

0. 0948 

0.0549 

0.3880 

0.2220 

0.  0898 

0.0566 

0.3980 

0. 1920 

0. 0818 

0.0479 

0.5040 

0.2240 

0. 0844 

0.0552 

0.3400 

0.2427 

0. 1054 

0.0649 

0.2300 

0.1907 

0.  0780 

0.0465 

0.3120 

0.1720 

0. 0802 

0.0506 

0.2820 

0.1793 

0.  0924 

0.0592 

0.2260 

0.2093 

0. 1070 

0.0622 

0.2260 

0.2093 

0.  1070 

0.0622 

0.2820 

0.1793 

0.  0924 

0.0592 

0.3400 

0.2427 

0. 1054 

0.0649 

0.2300 

0.1907 

0. 0730 

0.0465 
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TABLE    XXIX 


Normalized   Gradients   for    Flight    Track   Four 


5    meter 

15    meter 

50    meter 

100    mete 

0.5340 

0.2440 

0 .0948 

0.0549 

0  .  4  96  0 

0.2027 

0.0802 

0.0530 

0.3980 

0.1920 

0.0818 

0.0479 

0.5040 

0.2240 

0. 0844 

0.0552 

0.2300 

0.1907 

0.0780 

0.0465 

0.3880 

0.2220 

0.0898 

0.0566 

0.5040 

0.2240 

0  .0844 

0.0552 

0.2820 

0.1793 

0 .0924 

0. 0592 

0.3880 

0.2220 

0.0898 

0.0566 

0.2260 

0.2093 

0  .1070 

0.0622 

0.3400 

0.2427 

0.1054 

0.0649 

0.3120 

0. 1720 

0.0802 

0. 0506 

0.3120 

0.1720 

0.0802 

0.0506 

0.3400 

0.2427 

0.1054 

0.0649 

0.2260 

0.2093 

0 .1070 

0.0622 

0.3980 

0. 1920 

0.0818 

0. 0479 

0.2820 

0. 1793 

0.0924 

0.0592 

0.5340 

0.2440 

0 .0948 

0.0549 

0  .  4  96  0 

0.2027 

0.0802 

0.0530 

0.5340 

0.2440 

0.0948 

0.0549 

0.2300 

0.1907 

0.0780 

0. 0465 

TABLE    XXX 


Normalized   Gradients    for   Flight   Track    Five 


(Center) 


5    meter 

15    meter 

50     meter 

0.2840 

0.1253 

0.  0440 

0.1560 

0.0953 

0.  0448 

0.4560 

0.1807 

0.  0660 

0.3600 

0.  1507 

0.  0648 

0.3060 

0. 1773 

0.  0742 

0.1620 

0.1380 

0.  0602 

0.2840 

0.1627 

0.0840 

0.2820 

0. 1447 

0.  0670 

0.270  0 

0.1820 

0. 0974 

0.2100 

0.1593 

0.  0934 

0.0080 

0.0853 

0.  0716 

100  meter 

0.0249 
0.0243 
0.0382 
0.0337 
0.0436 
0.0414 
0.0520 
0. 0451 
0.0546 
0.0597 
0.0448 
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TABLE   XXXI 

Normalized  Gradients  for    Flight  Track    Five                (South) 

5    me+er           15   meter        50    meter  100    meter 

0.1707              3.06U6  0.0342 

0.1133              0. 0480  0.0262 

0.1253              0.0616  0.0355 

0.1687              0.0668  0.0461 

0.1733              0.0696  0.0489 

0.2407              0.0840  0.0476 

0.1707              0.0840  0.0512 

0.1400              0.0818  3.0522 


0 

.2460 

0, 

.2860 

0, 

.3380 

0. 

.2000 

0. 

.  1840 

0  , 

.3580 

0. 

.3240 

0. 

.2140 

TABLE    XXXII 

Normalized   Gradients  for   Flight    Track    Six                    (Center) 

5    meter          15    meter  50    meter        100    meter 

0.1220  0.0404  0.0232 

0.1507  0.0582  0.0283 

0.1713  0. 0694  0.0375 

0.1640  0.0684  0.0342 

0.1633  0.0734  0.0425 

0.1453  0.0680  0.0430 

0.2033  0. 0800  0.0483 

0.2033  3.0920  0.0544 

0.1400  0.0850  0.0588 

0.2113  0.0758  0.0440 

0.0780  3.0638  0.0421 


0, 

.2920 

0. 

.3200 

0. 

.3520 

0, 

.2220 

0  , 

.2780 

0. 

.2200 

0. 

.3540 

0, 

.3580 

0. 

.  1780 

0, 

.3540 

0. 

.  1680 
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TABLE    XXXIII 
Normalized  Gradients    for  Plight    Track    Six  (North) 

5    meter  15    meter        50    meter        100    meter 


0.2440 

0.106  0 

0. 0428 

0.0237 

0.2720 

0.1547 

0. 06  78 

0.0369 

0  .  4  06  0 

0.  1787 

0. 0730 

0.0377 

0.4480 

0.  1927 

0.  0758 

0.0423 

0. 1820 

0.1440 

0.  04  78 

0.0269 

0.3400 

0.1540 

0. 06  56 

0.0449 

0. 1900 

0.1313 

0.  0606 

0.0352 

0.1840 

0.1520 

0. 0686 

0.0489 

0.3280 

0.2280 

0. 0792 

0.0462 

0.3  72  0 

0.2007 

0. 0774 

0.0471 

0.2240 

0.0960 

0. 0456 

0.0359 
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APPENDIX  D 
MODEL  EQUATIONS  FOR  EACH  FLIGHT 

TABLE  XXXIV 
Model  Equations  for  MLD  (15  November) 

Model  Variables    Equations 

SST  MLD  = 

SST,    Position  MLD   =    -    1  22 

-    0.0 
+    14.6     (SST) 


-   0  .0918    {Position) 
"ST- 


SST.    Position.      MLD   =         -105 

Gradient    (100)  -   0.0718     (Position) 

+    10.5    (SST) 
♦    6  2.3    (GRAD(100) ) 

Not  es  * 

1.  SST    is    sea   surface  temperature. 

2.  Position   refers   to   station    position 
relative   to    the   most   seaward   station. 

3.  GEAD    (100)    refers    to   the   gradient 
computed    from   the    upper    100    meters    of 
the   thermocline. 
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TABLE   XXXV 
Model  Equations    for   MLD    (17    November) 

Model    Variables        Equations 

SST  MLD    =    +    43.7    +    1.01      (S5T) 

SST ,  Position     MLD  =  -  5  7.7  r      .  . 

-  0  .0865  (Position) 
+  10.9  (SST) 


SST  e  Position, 
Gradient  (15) 


MLD  =  ♦  62.8 

-  0  .098  (Position) 
♦  1 2.3  (SST) 

-  43.  1  (GRAE>(15)  ) 


SST.  Position. 

Gradient  (100)     MLD  =  -  8  8.2 


-  0  . 1 15  (Position) 

+  17.3  (SST) 

+  8  18  (GRAD  (100)  ) 


TABLE  XXXVI 
Model  Equations  for  MLD  (19  November) 

Model  Variables    Equations 

SST  MLD  =  +  68.7  -  0.86  (SST) 


SST,    Position  MLD    =    +•    17.9 

(SS"T) 


-   0.0551     (Position) 
+    3 .71     (SST) 


SST,    Position, 

Gradient    (5)  MLD    =    -   7.46 


-   0.045     (Position) 
'SST) 
GRAD(5)  ) 


+    4  .58     (SST) 
+   31.2    (r- 


SST  ,    Position. 

Gradient    (100)  MLD    =   ♦    18.2 


-  0.055     (Position) 
+    3  .72     (SST) 

-  8 .62    (GRAD(100) ) 
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TABLE  XXXVII 
Model  Equations  for  MLD  (1  December) 

Model  Variables   Equations 

SST  MLD  =  ♦  92.5  -2.03  (SST) 

SST,    Position  MLD    =    +    2  8.9 

-  0.061     (Position) 
+    4  .93     (SST) 

SST.  Position, 

Gradient  (50)      MLD  =  +  11.3 

-  0  .06  (Position) 
+  4  .62  (SST) 

♦  235.  (GRAD  (50) ) 

SST.  Position. 

Gradient  (100)     MLD  =  +  16.3 

-  0  .057  (Position) 
+  4  .28  (SST) 

+  344  (GRAD  (100)  ) 


TABLE  XXXVIII 
Model  Equations  for  MLD  (3  December,  Center) 

Model  Variables   Equations 

SST  MLD  =  +  9  5.5  -  2.48  (SST) 


SST,    Position  MLD    =   *   10.61 

-  0  .068 
*4.10  (SST) 


-  0  .0683  (Position) 
"ST* 


SST.  Position, 

Gradient  (15)      MLD  =  ♦  5  2.0 

-  0 .043 


-  0  .043  (Position) 

+  0 .95  (SST) 

+  1  18.2  (GRAD  (15)  ) 


SST.  Position. 

Gradient  (100)     MLD  =  +  7  1.7 


-  0  .064  (Position) 

-  3  .61  (SST) 

+  1384  (GRAD(100)) 
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TABLE    XXXIX 


Model   Equations   for    MLD    (3   December,    South) 


Model   Variables 

ssr 

SST,    Position 


SST  , .  Position, 
Gradient    (SO) 


SST,    Position. 
Gradient    (100) 


Equations 

MLD    =    +    95.5    -    2.48     (SST) 

MLD    =    ♦    40.61 

-  0  .0683     (Position) 

♦  4.10    (SST) 

MLD    =    +    1 27 

-  0.021     (Position) 

♦  0  .85     (SST) 

-  936.6     (GRAD  (50)  ) 


MLD    =    +  1 02 

-  0 .044  (Position) 

-  1  .530  (SST) 

♦  1  03.0  (GRAD  (100)  ) 


TABLE    XL 


Model   Equations   for    MLD    (5    December,    Center) 


Model   Variables 

SST 

SST,    Position 

SST,    Position, 
Gradient    (5) 


SST.    Position. 
Gradient    (100) 


Equations 

MLD    =    +    111    -    3.3  8     (SST) 

MLD    =    +    1 2.40 

-    0  .  103     (Position) 
+    8.18     (SST) 

MLD    =    +    8 .23 


-   0.10    (Position) 

+    7 .92    (SST) 

+    20.9     (GRAD(5)  ) 


MLD    ■    +  2  0.2 

-  0.099     (Position) 

♦  6.70     (SST) 

+  1  8"7.4     (GRAD  (100)  ) 
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TABLE   XLI 


Model  Equations   for    MLD     (5   Dscember,    North) 


Model   Variables 

SST 

SST,    Position 


SST,    Position, 
Gradient    (5) 


SST.    Position. 
Gradient    (100) 


Equations 

MLD    =+    1 23   -    5.01     (SST) 

MLD    =    -    4  3.7 

-0.12    (Position) 
+    13.1       SST) 


MLD    =    +    3.23 

-   0.10  (Position) 

♦    7.92  (SST) 

+    20.9  (GRAD(5)  ) 


MLD    =    -    34.0 

-   0  . 1 14     (Posi-ion) 
+    1  1.42      (SST) 

lD  (1 


+    1 74.7     (GRAD  (100) ) 
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